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I.  INTRODUCTION 


A.  BACKGROUND 

Surveying  is  the  science  and  art  of  measurements  which  are  necessary  to 
determine  the  relative  position  of  points  above,  on,  or  beneath  the  surface  of  the  earth, 
or  to  establish  the  points  in  a  specified  position.  Surveying  operations  are  conducted 
not  only  on  land,  but  also  in  the  oceans  and  in  space.  The  measurements  of  surveying 
consist  of  distances,  horizontal  and  vertical,  and  directions.  In  order  to  provide  a 
framework  of  survey  points  whose  horizontal  and  vertical  positions  are  accurately 
known,  basic  horizontal  and  vertical  control  surveys  must  be  performed.  A  primary 
use  of  control  surveys  is  for  construction  of  control  for  a  map  or  chart.  The 
fundamental  network  of  points  whose  horizontal  positions  have  been  accurately 
determined  is  called  horizontal  control  (Schmidt,  1978,  p.  122). 

Horizontal  control  generally  is  established  either  by  traverse,  triangulation,  or 
trilateration.  Which  one  is  to  be  used  depends  on  the  accuracy  required  and  the  factor 
of  economy  in  the  selection  of  survey  method.  Obviously,  there  are  many  degrees  of 
precision  possible  in  any  measurement  because  no  surveying  measurement  is  exact. 
Each  of  these  methods  may  be  the  best  one  to  use  for  a  given  purpose.  Ordinarily,  it  is 
a  waste  of  time  and  money  to  attain  unnecessarily  high  accuracy.  On  the  other  hand, 
if  the  measurements  are  not  sufficiently  precise,  faulty  survey  results  are  produced. 
Therefore,  the  best  surveyor  is  not  the  one  who  makes  the  most  precise  measurements, 
but  the  one  who  is  able  to  choose  and  apply  the  appropriate  measurement  with 
precision  requisite  to  the  purpose. 

Before  1950  the  main  framework  of  a  first-order  geodetic  survey  almost  always 
consisted  of  triangulation,  which  could  be  replaced  by  traverse  in  cases  where  the 
topography  made  triangulation  impracticable.  Today,  due  to  the  development  of 
electronic  distance  measuring  (EDM)  equipment,  the  first-order  control  points  can  be 
established  by  means  of  high  accuracy  traverse  [Allan,  1968,  p.  370].  Therefore,  the 
horizontal  control  is  frequently  provided  by  traverse,  especially  for  surveys  in  area  of 
limited  extent,  mostly  fiat  and  jungle  covered.  Traverse  in  such  cases  is  much  more 
economic,  convenient,  and  rapid  than  other  methods  and  the  results  are  equally 
accurate. 
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In  order  to  achieve  high  precision  of  horizontal  control  points  when  distributed 
over  a  large  area,  first  or  second-order  geodetic  surveys  are  required.  These  types  of 
survey  treat  the  shape  of  the  earth  as  ellipsoidal  and  require  using  the  most  accurate 
distance  and  angle  measurements.  Computation  of  such  a  survey  is  relatively 
complicated,  based  on  long  geodetic  formulas  for  computing  (with  necessary  precision) 
the  exact  horizontal  and  vertical  position  of  widely  distributed  points  on  the  earth's 
surface. 

Disregarding  ellipsoid  shape,  a  third-order  survey  is  used  over  earth  areas  of 
limited  extent.  In  this  type  of  a  survey,  the  earth  can  be  considered  as  flat  and  all 
angles  are  considered  to  be  plane  angles.  Surveys  of  this  type  are  used  in  the 
densification  of  geodetic  control. 

B.  OBJECTIVES 

As  mentioned  above,  the  traverse  method  has  been  used  worldwide  mostly  for 
densification  of  control  stations.  However,  there  are  many  methods  of  traverse 
computations.  The  main  objectives  of  this  thesis  are  to  (1)  compute  a 
closed-connecting  traverse  and  adjust  station  positions  by  the  Approximation  Method 
with  the  compass  (Bowditch)  rule  and  Least  Squares  Method  (adjustment  by 
observation  equations  only),  and  (2)  compare  the  results  of  the  two  methods. 

All  computations  have  been  accomplished  in  the  Universal  Transverse  Mercator 
(UTM)  grid  coordinates  rather  than  geodetic  coordinates. 
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II.  TRAVERSES 


A.  GENERAL 

1'hc  word  traverse  generally  means  to  pass  across.  But  in  surveying,  this  word 
means  the  measurement  in  a  specified  sequence  of  the  lengths  and  directions  of  lines 
between  points  on  the  earth  whose  position  may  be  know  or  unknown.  Traverse  is  the 
most  widely  employed  method  for  dcnsification  of  local  horizontal  control.  Linear 
measurements  are  made  either  by  direct  observation  using  a  tape  or  EDM  equipment, 
or  by  indirect  observation  using  tachcomctric  methods.  The  angular  measurements  arc 
made  with  theodolite  or  transit.  In  this  thesis,  the  only  traverse  operations  considered 
arc  for  angles  measured  by  theodolite  and  distances  measured  by  precise  EDM 
equipment  or  tapes. 

Two  kinds  of  traverse  exist  in  surveying,  an  open  and  a  closed  traverse. 

I.  Open  Traverse 

An  open  traverse  normally  originates  at  a  point  of  known  position  and  docs 
not  return  to  the  starting  point  nor  docs  it  terminate  on  another  point  of  known 
position  (Eigurc  2.1).  Open  traverses  should  generally  not  be  used  because  they  can 
not  be  checked  for  errors. 


Eigurc  2.1  Open  Traverse. 


2.  Closed  Traverse 

Closed  traverses  can  further  be  sub-divided  as,  a  closed-loop  and  a 
closed-connecting  traverse. 

A  closed-loop  traverse  is  one  that  originates  and  terminates  on  a  single  point 
of  known  position,  thus  forming  a  closed  polygon  (f  igure  2.2).  This  type  of  traverse 
provides  an  internal  check  on  angles  but  no  check  on  systematic  errors  in  distance. 
Also,  if  the  starting  azimuth  (between  stations  0  and  I  in  Figure  2.2  )  has  an  error,  it 
causes  error  in  orientation  of  the  entire  traverse.  A  closed-loop  traverse  generally 
should  not  be  used. 


Figure  2.2  Closed-loop  Traverse. 


A  closed-connecting  traverse  is  one  that  begins  and  ends  on  two  different 
points  whose  horizontal  positions  have  been  previously  determined  by  a  survey  of 
higher  or  equal  accuracy  (Figure  2.3).  This  type  of  traverse  is  preferable  to  all  others, 
since  computational  checks  arc  possible  to  detect  systematic  errors  in  both  distance 
and  direction. 

B.  ANGULAR  AND  DIRECTIONAL  MEASUREMENTS 

The  position  of  traverse  points  is  determined  by  the  direction  and  distance  from 
the  starting  point.  To  obtain  the  direction  by  means  of  azimuth,  the  horizontal  angle, 
or  plane  angle,  must  be  measured  in  the  field.  Also,  the  determination  of  vertical 
angles,  or  zenith  distances,  may  be  required  to  reduce  slope  distances  to  horizontal 
distances. 
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Figure  2.3  Closed-connecting  Traverse. 

Commonly,  both  horizontal  and  vertical  angle  measurements  arc  accomplished 
with  a  transit  or  theodolite.  The  theodolite  is  employed  especially  for  surveys  of  high 
precision.  Two  types  of  theodolite  arc  a  repeating  theodolite  and  direction  theodolite.  - 
A  repeating  theodolite  reads  directly  to  20"  or  T  and  may  be  estimated  to  one-tenth  of 
the  corresponding  direct  reading.  A  direction  theodolite  usually  reads  directly  to  1" 
and  may  be  estimated  to  0.1"  | Davis  ct  ah,  1981,  p.  2 1 5 J .  In  general,  a  direction 
theodolite  is  more  precise  than  a  repeating  theodolite  and  with  it,  plane  angles  arc 
computed  by  subtracting  one  direction  from  another. 

In  all  types  of  traverses,  the  horizontal  angles  can  be  measured  by  one  or  more 
of  the  following  listed  angle  measurement  methods. 

1.  Interior  Angle 

Interior  angle  is  the  angle  measured  within  a  closed  figure  at  the  intersection 
of  two  lines  (Figure  2.4). 

2.  Deflection  Angle 

Deflection  angle  is  the  angle  measured  from  the  extension  of  the  preceding  line 
to  the  succeeding  line  (Figure  2.5).  Such  angles  must  be  identified  as  right  or  left  to 
express  whether  the  angle  is  turned  to  the  right  or  to  the  left  from  the  prccccding  line. 

3.  Angle  To  The  Right 

Angle  to  the  right  is  the  clockwise  angle  measured  from  the  preceding  line  to 
the  succeeding  line  (Figure  2.6) 
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Figure  2.4  Measuring  of  Interior  Angles. 


Figure  2.5  Measuring  of  Deflection  Angles. 

C.  LINEAR  MEASUREMENT 

Direct  linear  measurements  may  he  obtained  in  traversing  by  pacing,  odometer 
reading,  tachcomctry  (stadia),  subtense  bar,  taping,  and  EDM.  Of  these  methods, 
taping  and  EDM  arc  most  commonly  used  by  surveyors.  However,  EDM  equipment 
has  a  clear  superiority  over  traditional  taping  for  lines  in  excess  of  about  250  meters. 

Distances  measured  using  EDM  equipment  arc  subject  to  personal,  instrumental, 
and  natural  errors.  Personal  errors  include  misreading,  improper  centering  of  the 
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Figure  2.6  Measuring  of  Angles  to  the  Right. 

instrument  over  the  stations,  failing  to  exactly  center  the  null  meter,  and  incorrectly 
measuring  meteorological  factors  and  instrument  heights.  Instrumental  errors, 
expressed  in  terms  of  the  accuracy  of  the  instrument  specified  by  the  manufacturer, 
contain  two  parts.  For  example,  if  the  accuracy  of  an  instrument  is  designated  as  ± 

(10  ppm  +  5  mm),  the  constant  error  part  is  +  5  mm,  which  is  independent  of  the 
distance,  and  the  value  of  the  proportional  part  is  10  ppm  (parts  per  million)  which  is  a 
function  of  the  distance  measured.  Constant  error  is  most  significant  for  short 
distances.  For  very  long  distances  the  constant  error  becomes  negligible,  but  the 
proportional  part  is  important.  Natural  errors  such  as  refraction  arc  cause  by 
changing  of  atmospheric  conditions  along  the  measured  line  between  the  end  stations. 

D.  ACCURACY 

In  survey  adjustment,  a  deviation  from  the  ‘true’  value  is  considered  as  an 
observational  error  and  the  standard  error  designates  the  measure  of  accuracy  of  the 
observation.  The  meaning  of  an  accuracy  is  then  the  degree  of  conformity  or  closeness 
of  a  measurement  to  true  value. 

The  quality  of  traverse  operations  is  dependent  upon  the  accuracy  of  angular  and 
linear  measurements;  thus,  in  checking  the  accuracy  of  traverse  two  quantities  arc 
considered,  the  angular  misciosurc  and  the  linear  misclosurc.  Although  the  positional 
closure  (relative  accuracy)  is  an  indication  of  the  overall  quality  of  the  traverse  and  is 
used  for  traverse  classification,  it  docs  not  yield  information  on  the  precision  of  point 
location  determined  in  a  traverse  [Davis  et  al.,  p.  332]. 
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The  inherent  weakness  in  a  traverse  is  that  the  deviation  of  each  measured  line  is 
determined  by  a  single  series  of  angular  observations,  further,  any  error  in  any  angle 
will  affect  not  only  the  adjoining  line  but  all  subsequent  lines  to  a  greater  or  lesser 
extent  according  to  their  lengths  [Allan,  1968,  p.  371], 

Angular  misclosure  is  expressed  by  standard  error  of  the  measured  angle  times 
the  square  root  of  the  number  of  measured  angles. 

Linear  misclosure  is  commonly  expressed  as  a  ratio  of  total  misclosure  to  total 
length  of  traverse. 

Finally,  some  of  the  most  significant  features  of  traverse  classification  by  the  U.S. 
Federal  Geodetic  Control  Committee  (1984)  are  shown  in  Table  I. 

E.  ADJUSTMENTS 

Adjustment  of  a  traverse  is  carried  out  to  ensure  consistency  within  the  known 
positions  of  the  originating  and  terminating  stations  and  to  remove  inconsistencies  in 
observed  angles  and  distances  to  compensate  for  random  errors.  For  a  more  precise 
extended  traverse,  adjustments  made  on  the  basis  of  least  squares  arc  preferred.  But  a 
traverse  of  limited  extent  can  be  adjusted  by  simple  approximation  methods. 

1.  Approximation  Methods 

There  are  four  methods  for  traverse  adjustment  by  approximation. 

a.  Arbitary  Method 

This  method  does  not  conform  to  a  fixed  rule.  Rather,  the  linear  error  of 
closure  is  distributed  arbitarily  according  to  arbitary  preference  of  the  surveyor. 

b.  Transit  Rule 

Transit  rule  is  better  for  adjustment  of  the  traverse  where  the  angles  are 
measured  with  greater  accuracy  than  distances,  and  is  valid  only  when  the  traverse  lines 
are  parallel  with  the  grid  system  used  for  the  traverse  computations.  Corrections  are 
made  by  the  following  rules:  the  correction  in  latitude  for  any  station  is  equal  to  the 
multiple  of  latitude  in  that  section  and  total  closure  in  latitude  divided  by  the  sum  of 
all  latitudes  in  traverse,  and  the  correction  in  departure  is  equal  to  the  multiple  of 
departure  in  that  section  and  total  closure  in  departure  divided  by  the  sum  of  all 
departures  in  the  traverse  [Davis  et  al.,  1981,  p.  323]. 

c.  Compass  or  Bowditch  Rule 

This  method  is  suitable  for  adjustment  of  the  traverse  where  the  angles  and 
distances  are  measured  with  equal  precision  and  uses  the  following  rules:  the  correction 
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TABLE  I 

TRAVERSE  CLASSIFICATION 


Order 

Class 

First 

Second 

I 

Second 

II 

Third 

I 

Third 

II 

Station  spacing  not 
less  than  (km) 

10 

4 

2 

0.  5 

0.  5 

Minimum  number  of 
network  control 
points 

4 

3 

2 

2 

2 

Theodolite  least 
count 

0.  2" 

1.  0" 

1.  0" 

1.  0" 

1.  0" 

Direction  number  of 
positions 

16 

8  or  12 

6  or  8 

4 

2 

Standard  deviation  of 
mean  not  to  exceed 

0.  4" 

0.  5" 

0.  8" 

1.  2" 

2.  0" 

Rejection  limit  from 
the  mean 

4  » 

5" 

5" 

5" 

5" 

Azimuth  closure  at 
azimuth  check  point 
( second  of  arc) 

1.  7VN  30VN 

4.  5VN  10.  OVN 

12VN 

Position  closure  0.  04,/K  0.  08Vk  0.  2>/k  0.  4,/K  0.  8Vk 
after  azimuth  or  or  or  or  or 

adjustment*  1: 100000  1: 50000  1: 20000  1: 10000  1: 5000 

N  =  route  distance  in  km 
K  =  number  of  segments 

*The  expression  containing  the  square  root  is  designed 
for  longer  lines  where  higher  proportional  accuracy  is 
required  and  the  results  are  in  meter. 

The  closure  (  e.  g.  ,  1:50,000)  is  relative  error  of 
closure. 


in  latitude  for  any  station  is  equal  to  the  multiple  of  the  length  in  that  section  and  total 
closure  in  latitude  divided  by  the  total  length  in  the  traverse,  and  the  correction  in 
departure  is  equal  to  the  multiple  of  the  length  in  that  section  and  total  closure  in 
departure  divided  by  the  the  total  length  in  the  traverse  [Schmidt,  1978,  p.  150]. 
d.  Crandall  Method 

Crandall  method  is  a  rather  complicated  procedure  which  is  more  rigorous 
than  either  the  compass  or  transit  rules  but  suitable  for  adjusting  traverses  where  the 
linear  measurements  contain  larger  random  errors  than  the  angular  measurement.  In 
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this  method,  the  angular  error  of  closure  is  first  distributed  in  equal  portions  to  all  of 
the  measured  angles,  then  linear  measurements  are  adjusted  by  using  a  weighted  least 
squares  procedure  [Brinker,  1977,  p.  22S], 

2.  Adjustment  by  Least  Squares  Method 

The  method  of  least  squares  adjustment  is  based  upon  the  theory  of 
probability;  it  simultaneously  adjusts  the  angular  and  linear  measurements  to  make  the 
sum  of  the  square  of  the  residuals  (error)  a  minimum  [Brinker,  1977,  p.  228],  This 
method  can  be  used  for  any  type  of  traverse.  Because  of  the  availibility  of  fast 
computing  devices  at  the  present  time,  the  least  squares  method  is  being  widely  used. 
Further,  the  least  squares  solution  has  the  advantage  that  it  determines,  quite 
objectively,  a  unique  solution  for  a  given  adjustment  problem  [Clark,  1973,  p.  121]. 

In  general,  adjustment  is  needed  whenever  there  are  redundant  observations 
(more  observations  than  are  necessary  to  solve  the  required  unknowns).  As  an 
example,  to  determine  the  angles  of  a  plane  triangle,  only  two  observed  angles  are 
required  because  the  third  angle  can  be  obtained  by  subtraction  from  180°.  When 
three  angles  are  observed,  the  sum  of  them  will  not  be  equal  180°  due  to  error  in 
measurements.  Therefore,  these  three  angles  should  be  adjusted  to  fit  the  functional 
model. 

The  redundancy  may  be  interpreted  to  mean  that  among  n  observations  there 
exist  r  conditions  or  functions  (n  >  r)  that  must  satisfy  the  model. 

Let  n  be  a  number  of  observations  and  nQ  the  minimum  number  of 
observations  to  find  the  uniquely  solution  in  the  model,  then  redundancy  or  degree  of 
freedom  in  the  statistic,  r,  is 

r  =  n  -  n0  (2.1) 

Consequently,  there  are  r  redundant  observations,  which  can  also  give  a 
solution.  To  detect  the  error  in  each  observation,  the  best  estimated  or  the  most 
probable  value  must  be  defined  because  the  true  value  is  not  known  exactly. 
Statistically,  the  best  estimated  value  of  a  group  of  repeated  observation  is  the  average 
(arithmetic  mean). 

Once  the  difference  between  observed  value  (XQ)  and  estimated  value  (X£)  is 
determined,  the  adjusted  value  (Xa)  is  obtained  through  a  least  squares  solution,  then 
the  residual  (v)  can  be  expressed  as 
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V 


(2.2) 


and 


Xa  =  Xe  +  dx 


(2.3) 


where  dx  is  the  correction  to  estimated  value  to  obtain  the  adjusted  value. 

The  least  squares  adjustment  method  is  based  upon  the  criterion  of  the  sum  of 
the  squares  of  the  observational  residuals  must  be  minimum. 

When  observations  are  considered  as  uncorrelatcd  and  of  equal  precision  (with 
identity  weight  matrix),  the  least  squares  condition  can  be  expressed  as 

<I>  =  +  =  Svi^  =  minimum  (2.4) 

or  in  matrix  form 

(j>  =  yTy  =  minimum  (2.5) 

where  V  is  the  vector  of  residuals. 

In  uncorrelatcd  observations  with  unequal  precision,  such  as  distances  and 
angles  [Mikhail,  1981,  p.  68],  the  Equations  2.4  and  2.5  become 

0  =  WjVj2  +  W2V22  +  ...  +  wnvn2=  £wivi2  =  minimum  (2.6) 

or  in  matrix  form 

(I>  =  V^WV  =  minimum  (2.7) 

where  wj  is  the  i1*1  clement  of  the  diagonal  weight  matrix  W  and  vj  is  the  residual 
associated  with  the  corresponding  i1^  observation. 

Generally,  the  relative  weights  are  inversely  proportional  to  variance,  thus  the 
weight  matrix  is  the  inverse  of  cofactor  matrix,  Q  (when  it  is  square  and  nonsingular) 
and  defined  as 


W  =  Q'1 


(2.8) 
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where  the  elements  of  cofactor  matrix  Q  are 

%  =  (2-9) 

and 

qjj  =  CTij/o-Q2  (2.10) 

where  <Tj~  the  variance  of  the  itF  observation,  <7-  is  the  covariance  between  the  itF  and 
jin  observations,  and  <7qZ  is  variance  of  unit  weight  [Mikhail,  1981,  p.  67]. 

For  the  case  of  uncorrelated  weight  observations,  the  cofactor  matrice  will  be 
diagonal  with  all  off-diagonal  elements  being  equal  to  zero,  thus  the  diagonal  elements 
of  weight  matrix  in  this  case  are 

wii  =  l'%  =  V^ii2  (2-11) 

Generally,  there  are  two  types  of  equation  in  least  squares  adjustment: 
condition  and  observation  equations.  Condition  equations  include  one  or  more 
observations  but  observation  equations  include  parameters  and  only  one  observation. 

The  condition  as  well  as  the  observation  equations  involved  in  an  adjustment 
problem  can  be  linear  or  nonlinear.  However,  least  squares  treatments  are  generally 
performed  with  linear  functions,  since  it  is  rather  difficult  and  often  impractical  to  solve 
nonlinear  models  [Mikhail,  1978,  p.  108].  Consequently,  whenever  the  equations  in  the 
model  are  originally  nonlinear,  they  have  to  be  linearized.  A  method  of  series 
expansion,  especially  Taylor's  series,  is  often  used  to  obtain  linear  equations.  Only  the 
zero  and  first-order  terms  are  used  and  all  other  higher-order  terms  are  neglected. 

Thus,  the  linearized  form  for  the  general  case  of  m  functions  of  n  variables  becomes 

F  =  F°  +  JmnAx  (2-12) 

where  is  the  zero-order  terms,  when  x  =  xP,  and  Jnm  is  a  Jacobian  matrix  of 
coefficients  of  first  order  of  n  variable  (Appendix  A). 

The  choice  between  the  observation  equations  (indirect  observation)  and 
condition  equations  (observation  only)  techniques  will  depend  mainly  on  the 
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mathematical  model  of  the  problem  to  be  solved.  However,  the  final  answers  are 
always  the  same. 

a.  Adjustment  by  Observation  Equations 

The  method  of  adjustment  by  observation  equations  is  performed  with  both 
observations  and  parameters.  Therefore,  the  number  of  equations  is  equal  to  the 
number  of  observations.  Using  the  example  at  the  beginning  of  this  section,  if  three 
measured  angles  and  their  residuals  in  a  plane  triangle  are  assummed  to  be  ct,  P,  y,  Vj, 
V2,  and  v3,  respectively,  and  the  adjusted  value  of  those  angles  is  xa  j ,  xa2>  an<^  xa3 
then  the  three  condition  equations  in  the  normal  form  (zero  at  the  right-hand  side)  are 
(a  +  vi)-xai  =  0 

(P  +  v2>  *  xa2  =  0 
(Y  +  v3)-xa3  =  0 


by  using  Equation  2.3  for  adjusted  values 

V1  '  dxl  =  xeU  a 
v2  '  dx2  =  xe2'  P 
v3  •  dx3  =  xe3*  y 

letting 


V1 

-1 

0 

0 

dx-. 

V  = 

v2 

v3 

/  B  = 

0 

0 

-1 

0 

0 

-1 

,  X  = 

dxi 

dx§ 

then,  expressed  in  matrix  notation  as 


x 


x 


el 

e2 

e3 


Y 


V31  +  B33X3j  =  F31  and  the  general  form  of  the  adjustment  of  observation 
equations  become 


Vnl  +  ®nu^ul 


=  F 


nl 


(2.13) 


where  V  is  an  n  by  1  vector  of  residuals;  B  is  an  n  by  u  matrix  of  coefficients;  X  is  an  u 
by  1  unknown  vector  in  which  u  is  the  number  of  parameters;  and  f  is  an  n  by  1  vector 
and  equals  P  -  XQ,  in  which  P  is  evaluated  using  the  estimated  value  X£. 
b.  Adjustment  by  Condition  Equations 

The  method  of  adjustment  of  condition  equations  only  has  no  parameters 
included  in  the  condition  equations.  Thus,  the  number  of  condition  equations  is  equal 
to  the  number  of  redundancies.  From  the  example  above,  n=  3,  nQ  =  2,  and  r=  1, 
which  is  the  number  of  condition  equations  being  set.  In  this  case,  because  the  sum  of 
interior  angles  must  equal  180°,  the  single  condition  equation  can  be  expressed  as 
a  +  P  +  Y  +  vj  +  V2  +  v3  =  180° 
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v  1 +  v2  +  v3  =  1 80°  -  a  -  p  -  y 


A  =  |  1,1,1  |  , 


V  = 


v] 

v- 


Then,  the  general  form  of  this  technique  is 


| 180-a-p-y | 


ArnVnr  = 


rl 


(2.14) 


When  the  conditions  are  originally  linear,  the  vector  F  is  usually  written  in  terms  of  the 
given  observations  as 


rl 


-  Pd  -  ArnXo  nl 


(2.15) 


where  A  is  the  coefficient  matrix  V.,  P  is  a  constant  term  (see  Section  Il.E.a),  XQ  is 
observed  values,  r  is  redundancy,  and  n  is  a  number  of  observation  [Mikhail,  1976,  p. 
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III.  DATA  ACQUISITION  AND  REDUCTION 


A.  DATA  ACQUISITION 

Taverse  data  used  in  this  thesis  were  obtained  from  field  work  accomplished  from 
25  September  thru  9  October  1972  by  CAPT  Glen  R.  Schaefer,  NOAA  Corps,  and  Mr. 
Jim  D.  Shea,  National  Ocean  Service  (NOS),  utilizing  traverse  methods  in  Pinellas 
County,  Florida.  Only  the  first  15  of  40  occupied  stations  and  three  intersection  points 
will  be  used  for  analysis  (Figure  3.1).  The  two  pairs  of  known  stations  for  this 
closed-connecting  traverse  are  shown  in  Table  II. 

The  known  stations  were  observed  by  the  US  Coast  and  Geodetic  Survey  (now 
NOS)  and  adjusted  by  the  National  Geodetic  Survey  (NGS).  Station  Turtle  2  is  of 
first-order  and  the  other  three  stations  are  of  third-order. 

The  horizontal  angles  measured  by  the  method  of  angles  to  the  right  (Table  III) 
were  turned  with  a  Wild  T-2  theodolite,  according  to  the  specifications  of  third-order 
class  I  traverse,  by  starting  at  station  Tomlinson  and  using  Egmont  Key  Lt.  House  for 
a  backsight.  The  traverse  was  closed  on  Turtle  2  with  a  check  azimuth  to  Madeira 
Beach  Tank. 

The  slope  distances  were  measured  in  the  field  with  a  Model  76  Geodimeter  in 
feet  and  corrected  for  temperature  and  pressure.  Distance  measurement  by  the  Model 
76  Geodimeter  are  reported  to  have  an  accuracy  in  the  temperature  range  of  -20°C  to 
+  50°C  of  ±(1  ppm  +  1  cm)  with  a  resolution  of  1  mm.  [Schmidt,  1978,  p.  116].  all 
observed  distances  were  converted  to  meters  and  reduced  to  horizontal  by  the 
procedure  given  later  in  this  chapter.  Finally,  geodetic  distances  are  reduced  to  grid 
distances  by  applying  the  scale  factor  correction  (Table  IV). 

B.  DATA  REDUCTION 

1.  Computation  of  Ellipsoidal  Distances 

For  the  requirement  of  high  precision  in  the  first-order  traverse,  the  measured 
distances  obtained  by  EDM  equipment  are  first  corrected  for  atmospheric  conditions 
and  then  reduced  to  the  ellipsoid  (Figure  3.2)  by  the  equations 

S  =  2  Ra  Sin'1(S0/2Ra)  (3.1) 


22 


f  6 

A, 


IP  3 


A  =  known  station 
O  =  traverse  station 
IP  =  intersection  point 


Nor  ro  Scale 


Figure  3.1  Traverse  Layout. 
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TABLE  II 

DATA  OF  KNOWN  STATIONS 

Station 

Station 

Grid  coordinate  ( m.  ) 

No. 

Name 

X/Easting 

Y/Northing 

0 

Egmont  Key  Lt.  House 

326214. 833 

3054014.  779 

1 

Tomlinson 

329400. 420 

3063485. 483 

15 

Turtle  2 

325348. 292 

3076179.  297 

16 

Madeira  Beach  Tank 

322698. 706 

3076362. 405 

s0  =  {  [m2-(h2-h,)]  /  (1  +  hj/RgXl  +  h2/Ra)  } ^2 


(3.2) 


1/Ra  =  {(Cos2a)/M}  +  ((Sin2a)/N}  (3.3) 


M  =  a(l-e2)/(l-e2Sin2<p)3/2  (3.4) 


N  =  a  /  (1  -  e2Sin2<p)1/2  (3.5) 

where  m  is  the  measured  distance,  S  is  the  ellipsoidal  distance,  Sq  is  chord  distance,  hj 
and  h2  are  the  heights  above  the  ellipsoid  at  each  end  of  the  line,  Ra  is  the  radius  of 
curvature  of  the  chord  distance  m  at  an  azimuth  a,  a  is  geodetic  azimuth  clockwise 
from  north,  M  is  radius  of  curvature  in  the  plane  of  the  meridian,  N  is  radius  of 
curvature  in  the  prime  vertical,  a  is  a  semimajor  axis  of  the  reference  ellipsoid,  e  is  its 
eccentricity,  and  <p  is  the  mean  latitude  of  that  line.  [Torge,  1980,  pp.  48,49,50,51,180] 
2.  Computation  of  Geodetic  Distances 

For  a  lower-order  traverse,  the  measured  distance  can  be  reduced  to  mean  sea 
level  (MSL)  or  geoid  only.  Because  the  error  of  only  1  ppm  in  length  results  from  an 
error  of  6  m  in  separation  between  MSL  and  the  spheroid  [Bomford,  1980,  pp.  42,  342, 
345]. 
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TABLE  III 

OBSERVED  HORIZONTAL  ANGLES 


A)  For  Traverse 


At  Stn.  Stn. 

Angles 

No. 

Name 

D 

M 

S 

0 

Egmont  Key  Lt.  House 

1 

Tomlinson 

105 

21 

25 

2 

TN-01 

243 

39 

18 

3 

TN-02 

168 

10 

15 

4 

TN-03 

59 

55 

56 

5 

Ruscue 

291 

28 

29 

6 

RE-01 

160 

43 

42 

7 

RE-02 

269 

55 

50 

8 

RE-03 

92 

43 

29 

9 

RE-04 

178 

22 

51 

10 

RE-05 

182 

31 

19 

11 

RE-06 

196 

54 

42 

12 

RE-08 

168 

51 

46 

13 

RE-09 

161 

06 

42 

14 

RE-10 

236 

58 

14 

15 

Turtle  2 

78 

37 

24 

16 

Madeira  Beach  Tank 

B)  For  Intersection  Points  ( IP) 


IP  #1  * 


At  Traverse 

Angles 

Back  Stn. 

Stn. 

D 

M 

S 

RE-01 

RE-02 

196 

23 

46 

RE-04 

RE-05 

184 

34 

22 

RE-05 

RE-06 

186 

29 

15 

RE-06 

RE-08 

34 

43 

31 

RE-08 

RE-09 

2 

41 

22 

*  * 

RE-03 

RE-04 

175 

35 

01 

RE-04 

RE-05 

97 

12 

06 

RE-05 

RE-06 

2 

32 

22 

k  k  k 

RE-08 

RE-09 

194 

01 

22 

RE-09 

RE-10 

232 

11 

34 

* 

•k  k 


St.  Petersburg  BCH  CO  Tank 
St.  Petersburg  BCH  St.  Johns  CH  Tower 
***  Bay  Pines  Veterans  Administration  Hosp. 


<7 

(  ±) 
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TABU:  IV 

MFASURFD  DISTANCES  AND  GRID  DIS 

I'ANCFS 

At  St n. 

Measured 

Distances 

Grid  Distances 

a 

To  Stn. 

ft 

m 

ft 

m 

(  ±m) 

1 

2300.  98 

701.  340 

2300.  91 

701.  318 

0.  017 

2 

2605.  41 

794.  131 

2605.  28 

794.  090 

0.  018 

3 

4452. 48 

1357. 119 

4452. 31 

1357.  068 

0.  024 

4 

709. 56 

216.  274 

709.  54 

216.  267 

0.  012 

5 

5050. 37 

1539. 356 

5050.  21 

1539. 307 

0.  025 

6 

6620. 54 

2017.  945 

6620. 33 

2017.  880 

0.  030 

7 

1655. 17 

504. 497 

1655.  12 

504. 481 

0.  015 

8 

1605. 85 

489. 464 

1605. 80 

489. 448 

0.  015 

9 

2114.  68 

644.  556 

2114.  61 

644.  535 

0.  016 

10 

2362. 13 

719. 979 

2362. 05 

719. 956 

0.  017 

11 

1360. 16 

414.  578 

1360.  10 

414.  561 

0.  014 

12 

5060.  98 

1542. 590 

5060. 41 

1542. 415 

0.  025 

13 

6172. 26 

1881. 309 

6171.  67 

1881.  128 

0.  029 

14 

15 

6664. 72 

2031.  411 

6664.  51 

2031.  346 

0.  030 

m 

= 

measured  distance 

So 

= 

chord  distance 

S 

S 

= 

arc  distance  on 

\ 

ellipsoid 

s0  y 

h, 

= 

elevation  above 

ellipsoid  at 

station  1 

\  A* 

'h 

elevation  above 

ellipsoid  at  station 

2 

= 

radius  earth  curvature 

along  measured  line 

0 

Figure  3.2  1  Ellipsoidal  Distance. 

However,  the  process  of  reduction  requires  three  steps:  (1)  correct  slope 
distances  to  horizontal  distances,  (2)  reduce  horizontal  distance  to  geodetic  distances, 
and  (3)  change  the  geodetic  distances  to  grid  distances. 

1  he  slope  distance  data  used  included  correction  for  atmospheric  conditions. 
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In  a  plane  survey,  such  as  a  traverse,  there  are  two  considerations  for  the 
reduction  of  slopojdistances  to  horizontal  distances,  a  short  slope  distance  and  a  long 
slope  distance. 

a.  Slope  Reduction  for  Short  Distances 

Short  slope  distances  (  <  2  ini  or  3.3  km)  measured  by  using  IiDM 
instruments  separate  from  the  theodolite,  can  be  reduced  to  horizontal  with  a  simple 
trigonometry  process  as 

D  =  m  cos  0  (3.6) 

where  D  is  the  horizontal  distance,  m  is  the  slope  distance,  and  0  is  the  vertical  angle 
(Figure  3.3). 


Figure  3.3  Slope  Reduction  for  Typical  Triangle. 

The  horizontal  distance  also  may  be  determined  by  using  the  difference  in 
elevation  between  the  two  ends  of  the  line.  The  horizontal  distance  is 

D  =  (n^-d2)1'2  (3.7) 


in  which  d  is  the  difference  in  elevation  between  the  two  end  points.  The  heights  of  the 
EDM  instrument  and  reflector  above  the  survey  mark  must  be  observed,  and  d 
becomes 
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d  =  (E1  +  i1)-(E2  +  o2) 


(3.8) 


where  Ej  and  E?  are  the  elevation  at  each  end  of  the  line  respectively,  I  j  is  the  height 
of  the  instrument,  and  02  is  the  height  of  the  reflector.  Then,  expanding  the  right  side 
of  equation  3.7  with  the  binomial  theorem  yields 

D  =  m-(d2/2S  +  d4/8S3  +  ...)  (3.9) 

The  quantity  in  the  parenthesis  is  designated  slope  correction.  For 
moderate  slope  the  first  term  is  usually  adequate.  When  the  slope  distances  and 
vertical  angles  are  obtained  by  separated  EDM  equipment  from  theodolite,  the 
correction  of  the  vertical  angle  must  be  determined.  The  corrected  vertical  angle  0y  is 

eT  =  e0+  A0  (3.10) 

where  0q  is  an  observed  vertical  angle  by  theodolite,  and  A0  is 

A0  =  H  Cos  0Q  /  S  Sin  1"  (3.11) 

and 


H  =  (Hr-Ht)-(HrHe)  (3.12) 

where  Hr  is  the  height  of  the  reflector,  Ht  is  the  height  of  the  target,  Hj  is  the  height  of 
the  EDM,  and  He  is  the  height  of  the  theodolite  [Davis  et  al.,  1981,  pp.  103-104]. 
Equations  3.11  and  3.12  are  not  needed  when  the  slope  distances  and  vertical  angles 
are  obtained  simultaneously  by  using  an  EDM  transmitter  built  into  a  theodolite, 
b.  Slope  Reduction  for  Long  Distances 

Slope  reduction  for  long  distances  (>  2  mi  or  3.3  km)  involves  using 
vertical  angles  affected  by  curvature  and  refraction.  By  assuming  a  mean  radius  for  the 
earth  of  3959  mi  or  6371  km,  then  the  curvature  correction  (C),  expressed  as  an  angle 
in  seconds,  is  4.935"  per  1000  ft  or  16.19"  per  km  and  the  horizontal  distance,  D,  is 

D  =  m  Sin(9O°-0-C)  /  Sin(90°  +  C)  (3.13) 
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for 


6  =  (Y  +  P)/2  (3.14) 

where  y  and  P  are  the  vertical  angles  at  each  end  of  the  measured  line  [Davis  et  al.,  pp. 
106-107], 

When  a  single  vertical  angle  (y)  is  observed,  0  is  the  corrected  vertical  angle 
for  combined  results  of  curvature  and  refraction  (C&R),  then,  0  is  y  +  (C&R)'\  The 
C&R  correction  is  4.244"  per  1000  ft  or  13.925"  per  1000  m.  The  correction  of  C&R 
will  be  positive  when  the  vertical  angle  is  an  elevation  angle  and  negative  in  the  case  of 
a  depression  angle  [Davis  et  al,  1981,  p,  108]. 

3.  Reduction  of  Horizontal  Distances  to  Geodetic  Distances 

The  horizontal  distance  at  same  elevation  above  the  geoid,  must  be  reduced  to 
a  geodetic  distance.  This  can  be  done  by  the  equation 

D'  =  (R)(D)/(R+E)  (3.15) 

where  D'  is  the  geodetic  distance,  R  is  the  mean  radius  of  the  earth's  surface  at  that 
section,  D  is  the  horizontal  distance  at  elevation  E  above  the  geoid  [Davis  et  al.,  1981, 
p.  107], 

C.  GRID  DISTANCES 

The  traverse  computation,  based  on  the  UTM  grid  coordinate  system,  requires 
the  reduction  of  geodetic  distances  to  the  plane  of  the  projection  by  applying  the 
projection  scale  factor  and  grid  scale  constant.  Scale  factor  can  be  obtained  from  a 
graph  or  from  a  rigorous  formula  [Department  of  the  Army  Technical  Manual,  1958, 
pp.  3,4,9,17];  i.e., 

k  =  kQ  [  1  +  (XVIII)  q2  +  0.00003  q4  ]  (3.16) 

where  k  is  the  scale  factor  at  scale  working  on  the  projection,  kQ  is  the  central  scale 
factor  which  is  an  arbitary  reduction  applied  to  all  geodetic  lengths  to  reduce  the 
maximum  scale  distortion  of  the  projection  (for  UTM,  kQ  =  0.9996),  and  values  for  q 
and  (XVIII)  are  obtained  by  the  formulas  which  shown  in  Table  V. 
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TABLE  V 

SPECIFICATION  OF  PARAMETERS 

XVIII 

= 

1  +  e'  Cos2<p  .  1  .  1012 

2  k7 

a 

V 

(  1  -  e2  Sin2(p ) 1//2 

e2 

= 

( ecentricity )2  =  (  a2  -  b2 )  /  a2 

e2' 

= 

e2  /  (1  -  e2) 

q 

= 

0. 000001  E' 

E 

— 

grid  easting  =  E'  +  500,000  when  point  is 
east  of  central  meridian,  500,000  -  E 
when  point  is  west  of  central  meridian 

V 

= 

radius  of  curvature  in  the  prime  vertical 

a 

= 

semi-major  axis  of  the  ellipsoid 

b 

— 

semi-minor  axis  of  the  ellipsoid 

<P 

— 

latitude 
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IV.  TRAVERSE  COMPUTATION  AND  ADJUSTMENT 


Linear  measurements  and  angles  must  be  checked  by  computation  to  determine 
the  position  of  traverse  stations  and  whether  the  traverse  meets  required  precision. 
Traverse  station  coordinates  are  usually  expressed  in  terms  of  geographical  coordinates 
or  rectangular  coordinates  such  as  those  based  on  an  UTM  projection.  Traverse 
computations  are  usually  done  in  rectangular  coordinates  because  of  the  ease  of 
computation.  In  this  thesis,  only  closed-connecting  traverse  computations  in  UTM 
coordinates  were  used.  When  specifications  were  satisfied,  the  traverse  was  adjusted 
for  perfect  geometric  consistency  among  angles  and  lengths. 

A.  DATA  PROCESSING 

1.  Set  up  of  data  base 

Two  files  on  IBM  370/3033AP  main  frame  computer  at  NPS  were  established. 

2.  Modification  of  an  Existing  Program 

The  Fortran  program  TRAVADJ,  originally  written  by  LCDR  Saman 
Aumchantr,  RTN  (1984)  in  Watfiv  language  for  computing  and  adjusting  the  traverse 
station  coordinates,  was  modified  to  be  able  to  handle  the  distances  reduction 
processes. 

3.  Writing  a  New  Program 

A  Fortran  program  INDTRA  was  written  for  computing  and  simultaneously 
adjusting  traverse  station  and  intersection  point  coordinates  by  least  squares 
observation  equations  method. 

B.  COMPUTATION  OF  STARTING  AND  CLOSING  AZIMUTHS 

The  directions  of  lines  by  means  of  azimuth  are  used  for  traverse  computation, 
because  sines  and  cosines  of  azimuth  angles  automatically  provide  correct  algebraic 
signs  for  latitudes  and  departures. 

The  terms  latitude  and  departure  are  widely  used  in  rectangular  coordinate 
calculations  of  surveying.  The  latitude  of  a  line  is  its  projection  on  the  reference 
meridian,  which  differs  from  geographic  latitude.  The  departure  of  a  line  is  its 
projection  on  the  east-west  line  perpendicular  to  the  reference  meridian.  In  traverse 
calculations,  north  latitudes  and  east  departures  are  considered  plus;  south  latitudes 
and  west  departure,  minus. 
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Latitudes  are  also  sometimes  termed  ‘northings’  and  ‘Y  dilFerences’  (AY); 
departures  are  similarly  called  ‘eastings’  and  ‘X  differences’  (AX). 

Because  the  closure  angle  of  traverse  can  be  checked  by  azimuth  of  each 
consecutive  line,  starting  and  closing  azimuths  have  to  be  determined  in  the  first  step 
of  computation  and  azimuths  can  be  computed  from  a  pair  of  known  coordinate 
station  positions  at  the  two  ends  of  the  traverse  (Figure  3.1) 

The  azimuth  of  the  line  from  A  to  B,  ctAg,  measured  clockwise  from  north,  is 
determined  by  the  equation 

aAB  =  Tan-*(AX/ AY)  (4.1) 


for 


AX  -  X„.  XA 


(4.2) 


and 


AY  -  YB-  ya  <4'3> 

where  XA  and  Xg  are  the  grid  easting  coordinates,  and  YA  and  Yg  are  the  grid 
northing  coordinates  of  stations  A  and  B,  respectively  (Figure  4.1).  The  quadrant  of 
the  azimuth  of  line  AB,  aAg,  is  dependent  upon  the  sign  of  AX  and  AY  (Table  VI). 
The  back  azimuth  dgA  (the  azimuth  from  B  to  A)  is  obtained  by  adding  180°  to  the 
forward  azimuth  «Ag. 

The  length  of  the  line  AB  (denoted  as  SAg  or  S)  can  be  determined  by  the 
Pythagorean  theorem  or  by  one  of  the  trigonometric  relationships 

S  =  AX  /  Sin  a  (4.4) 


or 


S  =  AY  /  Cos  a 


(4.5) 
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Figure  4.1  Azimuth  Computation. 


TABLE  VI 

QUADRANT  OF  AZ1MUT1 1 


Quadrant  of  azimuth  measured  Sign  of  Sign  of 
clockwise  from  north  AX  Af 


0  to  90  +  + 

90  to  180  + 

180  to  270 

270  to  360  -  + 


Substituting  data  from  Table  11  into  liquations  4.2  and  4.3,  the  AX  and  AY 
between  Fgmont  Key  Lt.  House  and  Tomlinson  arc  computed  as  AX  =  3185.585  m 
and  AY  =  9470.704  nr.  The  azimuth  of' the  line  from  Fgmont  Key  Lt.  House  to 
Tomlinson  is  18°  35'  27.6"  (Equation  4.1).  Similarly,  the  azimuth  Trom  Turtle  2  to 
Madeira  Beach  lank  is  computed  to  be  273°  57'  12.0".  These  starting  and  closing 
azimuths  will  be  used  Tor  computing  the  coordinates  of  each  traverse  station  and  lor 
checking  the  angular  error. 
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C.  COMPUTATION  OF  TRAVERSE  STATION  COORDINATES 

Computation,  of  traverse  station  coordinates  is  the  reverse  process  of  finding 
azimuth  and  distance  from  coordinates.  Therefore,  the  rectangular  coordinates  for 
each  closed-connecting  traverse  station  can  not  be  computed  unless  forward  azimuth 
and  distance  from  the  previous  station  are  known. 

The  azimuth  is  reckoned  clockwise  from  north  and  obtained  by 

ajk  =  a-  +  180°  +  Pj  (4.6) 

where  oijk  is  the  forward  azimuth  from  station  j  to  station  k,  a-  is  the  forward  azimuth 
from  station  i  to  station  j,  and  Pj  is  the  horizontal  angle  at  station  j  for  j  values  of  1  to 
n,  where  the  previous  station  i  =  j  -  1  and  the  next  station  k  =  j  +  1.  The  number  j 
will  increase  from  1  (which  designates  the  starting  known  station  of  the  traverse)  to 
number  n  (which  was  the  last  occupied  and  known  station  of  the  traverse). 

Departures  and  latitudes  are  then  computed  by  using  Equations  4.4  and  4.5 
which  are  rewritten  as 

(4.7) 

and 


^^jk  Sjk^n  ajk 


AYjk  -  SjkCos  «jk  (4.8) 

where  Sjk  is  the  distance  between  stations  j  and  k.  The  coordinates  of  all  other 
traverse  stations  can  be  determined  by  adding  successive  departures  (AXjk)  and 
latitudes  (AYjk)  to  the  previous  station's  X  and  Y  coordinates,  respectively. 

Using  the  data  in  Tables  II,  III,  and  IV,  the  azimuth  and  coordinate 
computations  at  the  first  station  are  shown  here  by  using  Equations  4.6,  4.7,  and  4.8. 
The  result  for  other  stations  can  be  seen  in  Table  VII. 

Computing  the  azimuth  from  station  1  to  2 

Az  =  18°35'27.6"  +  ISO0  +  105°21'25"  =  303°56'52.6" 

Computing  the  coordinates  at  station  2 

X-easting  =  329400.420  + [701.318  Sin(303°56'52.6")] 

=  328818.645  m 
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Y-northing  =  3063485.483  +  [701.3 18  Cos(303°56'52.6")] 

=  3063877.126  m 

When  the  coordinates  of  all  stations  have  been  computed,  they  are  still 
unadjusted  coordinates  and  can  then  be  adjusted  by  one  of  the  two  methods  of  Section 
II.E  (approximation  and/or  least  squares  method). 


TABLE  VII 

UNADJUSTED  TRAVERSE  STATION  POSITIONS 


Stn.  Angles  Forward  Dist.  Unadjusted 

Azimuths  (m)  Coordinates 

D  M  S  D  M  S  X(m)  Y(m) 


18 

35 

28 

1 

105 

21 

25 

*329400. 420 

*3063485.  483 

303 

56 

53 

701.  318 

2 

243 

39 

18 

328818. 645 

3063877.  126 

7 

36 

11 

794.  090 

3 

168 

10 

15 

328923. 709 

3064664.  236 

355 

46 

26 

1357.  068 

4 

59 

55 

56 

328823. 700 

3066017. 613 

235 

42 

22 

216.  267 

5 

291 

28 

29 

328645.  029 

3065895.  760 

347 

10 

51 

1539. 307 

6 

160 

43 

42 

328303. 493 

3067396.  699 

327 

54 

33 

2017.  880 

7 

269 

55 

50 

327231.  464 

3069106.  259 

57 

50 

23 

504.  481 

8 

92 

43 

29 

327658. 538 

3069374.  789 

330 

33 

52 

489. 448 

9 

178 

22 

51 

327418. 001 

3069801.  054 

328 

56 

43 

644. 535 

10 

182 

31 

19 

327085. 512 

3070353.210 

331 

28 

2 

719. 956 

11 

196 

54 

42 

326741. 616 

3070985.  723 

348 

22 

44 

414. 561 

12 

168 

51 

46 

326658. 106 

3071391.  785 

337 

14 

30 

1542. 415 

13 

161 

6 

42 

326061.  428 

3072814.  113 

318 

21 

12 

1881.  128 

14 

236 

58 

14 

324811.  349 

3074219.  797 

15 

19 

26 

2031.  346 

15 

78 

37 

24 

273 

56 

50 

325348. 180 

3076178.  924 

*  Coordinates  for  station  1  were  known  and  held  fixed. 
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D.  ADJUSTMENT  BY  APPROXIMATION  METHOD 

In  this  thesisT  the  method  of  Compass  or  Bowditch  rule  was  used  to  adjusted  the 
data  in  Tables  III  and  IV.  Thus,  the  first  step  is  to  determine  the  angular  error  of 
closure  and  adjust  the  angles  to  obtain  the  proper  closing  azimuth  (closed  azimuth  at 
fixed  stations). 

1.  Angular  Errors  of  Closure 

In  the  closed-connecting  traverse,  when  there  are  n  stations  of  observed 
horizontal  angles,  n-1  lines  will  be  measured.  An  angular  error  in  traverse  can  be 
checked  and  obtained  at  the  last  station  by  comparing  the  computed  azimuth  and 
closing  azimuth  at  the  known  station. 

The  closing  azimuth  computed  (from  the  known  station  coordinates  at  the 
traverse  end)  at  the  station  15  is  273°  57'  12.0".  But  because  of  error  in  measurement, 
the  azimuth  computed  through  the  traverse  at  this  station  is  273°  56'  49.6",  which  is  a 
difficiency  of  22.4".  This  amount  of  angular  error  in  15  observed  stations  meets  the 
limit  for  allowable  error  for  a  third-order  class  I  traverse  (allowable  error  from  Table  I 
is  38.7"). 

The  average  correction  (Table  VIII)  is  distributed  uniformly  over  all  the  15 
traverse  angles  (Table  III). 

2.  Linear  Errors  of  Closure 

When  all  angles  have  been  corrected,  the  process  of  calculating  the  improved 
coordinates  of  all  traverse  stations  may  be  done.  The  check  on  angular  closure  for  a 
closed  traverse  does  not  guarantee  that  the  entire  survey  is  correct,  because  there  can 
be  considerable  errors  in  the  linear  measurement  of  individual  lines.  Such  errors  may 
not  show  up  in  the  angular  check.  In  order  to  check  the  closure  of  the  traverse,  it  is 
necessary  to  determine  linear  error. 

The  linear  error  (LE),  the  departure  error  (5x),  and  latitude  error  (Sy)  in  a 
traverse  are  determined  by  equations 

LE  =  [  (Sx)2  +  (5y)2  ]E 2  (4.9) 


5x  = 


GEn  -  GEn 


(4.10) 
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(4.11) 


6y  -  GNn  ■  GNn 


where  5x  and  Sy  are  the  traverse  closure  in  departure  and  latitude,  GEn  and  GEn'  are 
the  known  and  computed  grid  easting,  and  GNn  and  GNn'  are  the  known  and 
computed  grid  northing  at  the  closing  station,  respectively.  By  substituting  the  fixed 
and  computed  values  of  grid  easting  and  northing  in  Equations  4.10  and  4.11  for  the 
data  of  Tables  II  and  VII 

6x  =  325348.292  -  325348.180  =  +  0.112  m 

Sy  =  3076179.297  -  3076178.924  =  +  0.373  m 
LE  =  [(0. 1 1 2)2  +  (0.373)2] 1/2  =  0.389  m 

The  relative  error  of  closure  provides  a  better  assessment  of  the  quality  of  a 
traverse  than  the  linear  error  of  closure.  Therefore,  it  is  common  practice  to  calculate 
the  relative  error  of  closure,  which  is  the  linear  error  of  closure  divided  by  total 
distances  of  traverse,  and  to  express  the  result  in  the  form  of  a  ratio  with  unity  as  the 
numerator.  For  the  data  of  Tables  II  and  VII,  this  computation  is  0.389  /  14853.800 
or  1  :  38,185. 

Using  the  Compass  or  Bowditch  rule,  the  computed  traverse  closures  and 
corrections  were  obtained  (Tables  VIII  and  IX)  and  then  the  adjusted  station 
coordinates  (Table  X)  were  computed. 


E.  LEAST  SQUARES  ADJUSTMENT  BY  OBSERVATION  EQUATIONS 

The  adjustment  by  observation  equations,  shown  in  general  form  by  Equation 
2.13,  can  be  done  more  directly  than  the  adjustment  by  condition  equations.  To 
achieve  this,  the  explicit  expression  for  the  residual  vector  V  from  Equation  2.13  is 
substituted  in  Equation  2.7  to  obtain  the  following  equation: 

O  =  (F  -  BX)T  W  (F  -  BX)  (4.12) 

=  (FT  -  BTXT)(WF  -  WBX) 

=  XTBTWBX  -  XTBTWF  +  FT\VF  -  FT\VBX 
=  XTBTWBX  -  2FT\VBX  +  FTWF 

where  X^B^WF  =  F^WBX  are  scalar  quantities. 
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TABLE  VIII 
TRAVERSE  CLOSURE 


I)  Angular  error 


Known  azimuth  at  last  station 

Computed  azimuth 
Angular  error 

Angular  correction  per  station 


273°  57'  11.95" 

273°  56'  49.58" 
22. 37" 


II)  Linear  error  before  adjusting  azimuths 


Known  Coordinate 
at  Turtle  2 


X(m) 

325348.  292 


Y(m) 

3076179.  297 


Computed  Coordinate 

at  Turtle  2  325348.180  3076178.924 

Error  5x  =  +  0.  112  Sy  =  +  0.  373 


Linear  error  of  closure 
Total  distances 
Relative  error  of  closure 


0.  389  m 
14853.  800  m 
1  /  38, 185 


TABLE  IX 

LATITUDE  AND  DEPARTURE  CORRECTIONS 


Traverse  lines 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


X-correction  (m) 

-  0.  031 

-  0.  035 

-  0.061 
-  0.  010 

-  0.  069 

-  0.  090 

-  0. 022 
-  0. 022 

-  0.  029 

-  0.  032 

-  0.  018 

-  0.  069 

-  0.  084 

-  0.  090 


Y-correction  (m) 

+  0.  007 
+  0.  008 
+  0. 012 
+  0. 002 
+  0.  015 
+  0. 019 
+  0.  005 
+  0. 005 
+  0.  006 
+  0.  007 
+  0.  004 
+  0.  015 
+  0.  018 
+  0.  020 
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ADJUSTED  COORDINATES  BY  COMPASS  RULE 
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Coordinates  of  stations  1  and  15  were  known  and  held  fixed. 


The  minimization  of  Equation  4.12  can  be  done  by  taking  partial  derivative  with 
respected  to  each  ilnknown  variable  (X) 

=  d<t>/dx  =  2XTBTWB  -  2FTWB  =  0  (4.13) 

Transposing  Equation  4.13  and  rearranging  yields 

(BTWB)X  =  BTWF  (4.14) 


or 


NX  +  U  =  0  (4.15) 

T  T 

where  N  =  (B  WB)  is  the  normal  matrix  of  dimension  u  by  u,  and  U  is  B  WF.  Then 

X  =  -N_1U  (4.16) 

For  the  adjustment  in  traverse,  the  vector  V  in  equation  2.13  is  represented  the 
residual  of  observed  angles  and  distances.  If  there  are  n  observed  angles  in  a  traverse, 
there  will  be  n  -  1  observed  distances  and  the  number  of  residuals  becomes  2n+  1 
which  includes  the  residual  of  angles  and  distances. 

There  are  two  types  of  condition  equations  in  the  adjustment  of  observation 
equations:  the  angle  condition  and  distance  condition. 

From  Figure  4.2  the  angle  condition  can  be  expressed  by 


via  P;  *  (®ij  *  ®ik)  (4-17) 

where  vja  is  a  residual  of  observed  angle,  ay  is  a  forward  azimuth,  a^,  is  a  backward 
azimuth  and  Pj  is  an  observed  angle.  Equation  4.17  is  suitable  when  ay  ^  a^.  If  ay 
<  ajj.,  the  quantities  in  parenthesis  must  be  added  by  360°. 

The  Equation  4.17  can  be  expressed  in  coordinates  of  the  two  stations  as 


„  1  X.  -  X.  ,  X.  -  X. 

via  =  Pi  -  (  tan*  ^ TT )  -  tan^M— J )  ] 

j  '  1  i  1  k  *  *  i 


(4.18) 
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Figure  4.2  Determination  of  Angle  and  Distance  from  Coordinates. 


and 


Si  = 


[  tan'*  (— 1 — — M  -  tan**(— - -  )  | 

Y.  -  Y.  Y,  -  Y. 


X.  -  X. 


-1/VX- 


(4.19) 


J  •  k 

When  linearized  by  using  liquation  2.12,  Liquation  4.19  yields 


gj  =  gjtXj',  Yj',  Xj',  Yj',  Xk\  Yk')  +  a,5Xj  +  a^Yj  +  a35Xj  +  a4SYj  + 
a55Xk  +  a6«Yk  ‘  (4.20) 

where  Xj',  Yj',  Xj',  Yj',  Xk',  and  Yk'  arc  the  estimated  station  coordinates.  By 
substituting  the  linearization  of  function  gj,  liquation  4.18  can  be  written  as 


via  +  al6Xi  +  a25Yi  +  a35Xj  +  a45Yj  +  a56Xk  +  a65Yk  =  F 


(4.21) 


where  5Xj,  5Yj . oXk,  and  SYk  arc  unknown  parameters  (correction  in  X  and  Y 

coordinates),  cocfiiccints  aj,  a 2 . a^,  and  a^  arc  the  partial  derivative  of  function  gj 

with  respected  to  Xj.  Yj . .  Xk,  and  Yk,  respectively. 

For  the  distance  condition,  it  can  be  expressed  as 


VW  -  j;  •  I  (Xj  -  Xj)2  4-  (Yj.Yj)2!1''2 


(4.22) 
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where 


hi  =  [  (Xj  -  Xi)2  +  (Yj  -  Yj  )2  I1/2 
The  linerized  form  of  Equation  4.23  is  then  given  as 


(4.23) 


Xk\  Yk')  +  b,8Xi  +  b28Yi  +  b38X-  +  b48Y-  (4.24) 


Substituting  the  linearized  form  from  Equation  4.24,  Equation  4.22  becomes 


vid  +  bjSXj  +  b28Yi  +  b38X-  +  b48Y-  =  F2  (4.25) 

where  FI  and  F2  represent  the  constant  form  for  Equations  4.21  and  4.25.  Table  XI 
shows  the  partial  derivatives  for  Equations  4.21  and  4.25. 

Sometimes,  the  station  positions  which  is  determined  by  intersection  from 
traverse  station  is  also  to  be  adjusted  simultaneously  with  traverse  station  positions. 

In  this  case  only  the  angle  conditions  are  added  to  observation  equations  of  traverse. 

In  the  data  adjusted  under  this  thesis,  there  are  15  observed  angles  and  14 
oberved  distances  which  consisted  of  13  unknown  points,  consequently,  there  will  be  29 
observation  equations  included  26  parameters  of  8X  and  5Y.  Thus,  they  can  be 
expressed  in  matrix  form  of  Equation  2.13  as 


V29,l  +  B29,26  X26,l  “  F29,l 


(4.26) 


where  Vj  j,  v2  j, . v^  j  are  the  residuals  of  angles;  v^  j,  Vj  -j  j, . v^  j  are  the 

residuals  of  distances;  B  is  the  coefficeint  matrix  (consistings  of  a's  and  b's)  of 
parameter  X  (Table  XI);  and  F  is  the  constant  vector. 

When  three  intersection  points  with  10  observed  angles  were  added  (Figure  3.1) 
to  adjust  the  coordinates,  the  observation  equations  have  39  equations  including  32 
unknown  parameters. 


V39,l  +  B39,32  X32,l  “  F39,l 


(4.27) 
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TABLE  XI 

THE  COEFFICIENTS  OF  ANGLE  AND  DISTANCE  CONDITIONS 


a 

a 

a 

a 

a 

a 

b 

b 

b 

b 


1 

2 

3 

4 

5 

6 

1 

2 

3 

4 


dX± 

d*i 

5Yi 

dgj  

d<3i 

dYi  = 
dgl 
5xk  ” 

SYk 

dX± 

5h 


dh.^ 

W1 


=  + 


=  + 


=  + 


Yi  -  Yi  +  Yk  -  y> 
(Sij)2  (Sik)2 
X,  -  X,  Xlr  -  Xi, 


=  + 


=  + 


<sij)z 

Yi  -  Yi 

(Sij)2 
xi  -  xi 

(Sij)* 

Yk  -  Yi 

(Sik>2 
xk  -  xi 

(Sik)2 
xi  -  xi 

siJ 

Yi  -  Yi 

Sij 

Xi  -  Xi 

sij 

Yi  -  Yi 

Sij 


(sik> 


By  using  Equation  4.16,  N'*  is  the  cofactor  matrix  and  the  diagonal  terms  of  this 
matrix  gives  the  variances  of  the  adjusted  coordinates.  To  obtain  the  residuals  of  all 
observed  quantities,  the  reverse  process  must  be  done.  The  correction  of  X  and  Y 
adjusted  coordinates  and  their  standard  deviation  (<t's)  are  shown  in  Table  XII.  And 
the  adjusted  standard  deviation  of  observed  quantities  were  obtained  by  multiplying 
standard  deviation  of  unit  weight  [(V  ^WV  /  r)*^]  to  a  squares  root  of  diagonal 
element  of  B(BTWB)'1  BT  matrix  (Table  XIII). 
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TABLE  XII 

ADJUSTED  COORDINATES  BY  OBSERVATION  EQUATION 
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TABLE  XIII 

ADJUSTED  STANDARD  DEVIATION  OF  ANGLES  AND  DISTANCES 


Number 


Angles 

Di stances 

seconds 

meter 

(  ± ) 

(  ±  J 

9.  8 

0.  029 

8.  9 

0.  035 

9.  1 

0.  045 

9.  2 

0.  024 

9.  3 

0.  048 

9.  4 

0.  056 

8.  4 

0.  030 

8.  5 

0.  029 

7.  8 

0.  032 

7.  1 

0.  034 

9.  6 

0.  028 

9.  7 

0.  048 

9.  8 

0.  054 

6.  8 

0.  058 

6.  8 
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V.  ANALYSIS  OF  RESULTS 


A.  COMPARISON  OF  ADJUSTED  COORDINATES 

If  the  given  coordinates  of  the  control  points  are  assumed  to  be  error  free,  then 
the  accuracy  of  the  traverse  station  coordinates  depends  only  on  the  accuracy  of 
distance  and  angle  measurements.  The  adjusted  traverse  coordinates  obtained  by  the 
approximation  method  are  of  a  lower  order  of  accuracy  as  only  the  errors  in  the 
misclosure  in  azimuths  and  distances  were  determined.  These  errors  were  distributed 
by  assuming  that  all  observed  quantities  had  an  equal  probable  occurrence. 

The  least  squares  adjustment  method  provides  a  better  approximation  of  the  true 
value.  Therefore,  the  adjusted  traverse  coordinates  obtained  by  this  technique  provided 
better  estimates  for  position  of  all  traverse  stations  and  the  accuracy  of  the  adjustment 
can  be  checked  and  statistically  tested. 

After  the  13  adjusted  traverse  station  (from  stations  2  to  14)  coordinates 
obtained  by  the  approximation  method  and  by  the  least  squares  method  were 
compared,  the  difference  in  coordinates  at  each  station  were  computed  and  plotted 
(Table  XIV  and  Figure  5.1).  The  largest  difference  was  at  station  14.  Because  stations 
1  and  15  are  held  fixed,  the  least  squares  techniques  adjusts  simultaneously  errors  in 
azimuths  and  distances,  while  the  adjustment  by  approximation  adjusts  errors  in 
azimuths  and  distances  sequentially.  Consequently,  the  largest  difference  in  traverse 
distances  occurs  at  the  last  station  (station  14)  before  closing  of  traverse  at  the  fixed 
station. 

When  the  standard  deviation  of  observed  quantities  before  the  adjustment 
(Tables  III  and  IV)  were  compared  with  those  obtained  through  adjustment  (Table 
XIII),  the  standard  deviation  of  all  observed  quantities  in  Table  XIII  showed 
increments.  That  means,  the  estimated  standard  deviations  were  optimistic. 

In  this  thesis,  the  three  intersection  points  were  also  adjusted.  The  adjusted 
coordinates  of  these  were  compared  to  NOS  results.  The  largest  difference  occurs  at 
point  no.  3  (5x  =  +  0.140  m;  5y  =  -  0.158  m).  The  standard  deviation  of  adjusted 
coordinates  at  this  point  are  crxj  =  ±  1.87,  cry3  =  ±  0.96  m  (Table  XV). 
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TABLE  XIV 

COMPARISON  OF  ADJUSTED  COORDINATES/DISTANCES 
OBTAINED  BY  APPROXIMATION  AND  LEAST  SQUARES  METHODS 


Differences* 

Coordinates  Distances 


stn. 

dx 

2 

f>%9 

_ 

3 

- 

0.  015 

- 

4 

+ 

0.  004 

- 

5 

- 

0.  013 

- 

6 

+ 

0.  002 

- 

7 

+ 

0.  008 

+ 

8 

+ 

0.  001 

- 

9 

+ 

0.  005 

- 

10 

- 

0.  008 

- 

11 

- 

0.  003 

- 

12 

- 

0.  014 

+ 

13 

- 

0.  005 

+ 

14 

- 

0.  019 

- 

*  Approximation  minus 


dy 

(m) 

( m) 

o;  007 

0.  020 

0.  017 

0.  023 

0.  022 

0.  022 

0.  012 

0.  018 

0.  008 

0.  008 

0.  024 

0.  025 

0.  005 

0.  005 

0.  013 

0.  014 

0.  023 

0.  024 

0.  034 

0.  034 

0.  003 

0.  014 

0.  041 

0.  041 

0.  101 

0.  103 

least 

squares  solution 

B.  ANALYSIS  OF  THE  REFERENCE  VARIANCE  OF  UNIT  WEIGHT 

The  weight  matrix  was  set  for  the  least  squares  adjustment  by  using  Equation 
2.1 1.  The  <Tq  (a  priori  reference  variance  of  unit  weight)  was  assumed  as  1.  The 
result  of  Gq  =  (V  1  WV  /  r)  was  obtained  after  adjustment.  The  value  of  Cq  (a 
posteriori  reference  variance  of  unit  weight)  can  be  used  to  evaluate  the  weighting 
scheme  used  in  the  least  squares  adjustment. 

The  standard  deviation  of  the  observed  angles  (<ra)  and  distances  (<r^)  used  in  the 
least  squares  adjustments  (solutions  1,  2,  and  3)  and  the  corresponding  a  posteriori 
obtained  for  these  solutions  are  listed  in  Table  XVI.  As  the  a  posteriori  G()-  for  the 
solution  no.  3  is  closest  to  the  assumed  a  priori  <7q  (=  1),  the  weight  (or  the  standard 
deviations)  used  for  observed  angles  and  distances  in  this  case  seem  the  most  realistic. 
Further  statistical  testing  for  Gq  done  by  Chi-squares  or  F-test  was  not  earned  out 
under  this  thesis. 
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DISTANCE  DIFFERENCES  (M) 

0.00  0.02  0.04  0.06  0.08  0.10 


Oi 

* 


Figure  5.1  Comparison  of  Adjusted  Distances  Obtained 
by  Approximation  and  Least  Squares  Methods. 
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TABLE  XV 

COMPARISON  OF  COORDINATES  AT  INTERSECTION  POINTS 


X(m) 

crx(  ±m) 

Y(m) 

cry(  ±m) 

At 

Point  No. 
NOS 

Adjusted 

1 

326616. 356 
326616. 470 

1.  72 

3071294. 683 
3071294.  635 

1.  03 

8x 

-  0.  114 

Sy 

+  0.  048 

At 

Point  No. 
NOS 

Adjusted 

2 

327056. 079 
327056. 183 

1.  68 

3070340. 474 
3070340. 444 

1.  14 

8x 

-  0.  104 

5y 

+  0.  030 

At 

Point  No. 
NOS 

Adjusted 

3 

325378. 101 
325377. 961 

1.  87 

3077263.378 
3077263. 536 

0.  96 

5x 

+  0. 140 

5y 

-  0. 158 

TABLE  XVI 

COM  PARI  SION  OF  VARIANCES  OF  UNIT  WEIGHT 

Solutions 

Variance  used  in  Adjustment 

VTWV/r 

2 

/N  2 

ffa  ad  ff0 

*0 

1 

2"  (  lOppm  +  0.  5  cm)  1 

15.  08 

2 

5"  (lOppm  +  1. 0  cm)  1 

3.  98 

3 

10"  (  lOppm  +  2.0  cm)  1 

1.  32 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

By  using  the  weight  of  observed  quantities  for  the  adjustment,  the  traverse 
station  coordinates  computed  and  adjusted  by  the  least  squares  observation  equations 
method  were  more  accurate  than  those  obtained  by  the  approximation  methods.  Even 
though  the  observation  equations  method  may  require  a  greater  number  of  equations 
than  the  condition  equations  method,  the  processing  of  the  data  for  adjustment  is 
easier  and  the  corrections  in  X  and  Y  coordinates  can  be  directly  obtained  through 
iterative  solution.  This  method  is  suitable  when  a  computer  with  a  memory  capacity  of 
over  500  K  bytes  is  available.  However,  for  local  work  or  a  relatively  short  traverse,  an 
approximation  method  is  commonly  utilized  when  economic  and  logistic  criteria  are 
considered. 

B.  RECOMMENDATION 

The  INDTRA  Fortran  program  written  for  this  thesis  is  automated  for  handling 
only  two  kinds  of  survey  techniques:  traversing  and  intersection.  With  the  computers 
available  at  NPS,  the  development  of  adjustment  programs  for  covering  a  wide  range 
of  survey  techniques  should  be  done  to  use  and  continue  analysis  of  the  mixed  kind  of 
survey  techniques  including  traverse,  triangulation,  trilateration,  resection,  and 
intersection. 
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APPENDIX  A 
LINEARIZATIONS 

This  is  a  linearized  form  of  m  functions  in  n  unknowns. 


y  l  =  r1(x1,x2,...,xn) 

y2  =  f2(x1,x2,..,xn) 


ym  =  Uxl’x2’-xn) 


yO  = 


y.° 

f^Xj0,  x2°, ... 

-  xn°) 

>2° 

f,(x2°,  X2°,  .. 

-  Xn°) 

ym° 

f  (x  X  °, 

mv  m  ’  m  ’ 

Xn°> 

dY 

~dX 


8_h 

dh 

<5y, 

^Xj 

dx2 

8xn 

dy 

J  m 

dy2 

5xn 

AX  = 


Ax. 

AxJ 


The  general  form  of  linearized  functions  becomes 
Y  =  Y°  +  JyxAX 
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APPENDIX  B 

TRAVADJ FORTRAN  PROGRAM 


This  program  is  used  for  computing  and  adjusting  the  traverse  station  position  by 
approximation  method  (Compass  rule). 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


FORTRAN  PROGRAM  "TRAVADJ" 

THIS  PROGRAM  IS  USED  FOR 

1)  REDUCE  SLOPE  DISTANCE  TO  ELLIPSOIDAL  DISTANCE 

2)  DETERMINE  GRID  DISTANCE 

3)  COMPUTE  CLOSE  TRAVERSE 

4)  ADJUST  COORDINATE  BY  COMPASS  RULE 


n  a  ti 

2)  VALUE  OF  1/F  (EX.  1/F  =  294.978698 

3)  CENTRAL  SCALE  FACTOR 

4)  LAT°  AT  STARTING  AND  CLOSING  POINT 

"idl , iml , si ,  id2 . im2 , s2" 

5)  GRID  N.  AND  L  OP  4  fcNOWN  STATION 

"qnl ,gel ,qn2,qe2,qn3 ,ge3.gn4,gn4" 

6)  NUMBER  OF  MEASURE  DIStANC£ 

"n 

7)  ELEVATION  AT  FIRST  OCCUPIED  POINT 

8' 

9 

10 
11 


12) 


)  NAME  OF  ALL  STATIONS 
)  INDICATOR  VALUE  1  =  NO  VERTICAL  ANGLE 
'  2  =  VERTICAL  ANGLE 

)  DIFFERENT  IN  ELEVATION  BETWEEN  TWO  STATIONS 
)  SLOPE  DISTANCE  (WITH  UNIT  FEET  OR  METER) 

HORIZONTAL  ANGLES 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

INPUT  DATA 

1)  SEMIMAJOR  AXIS  OF  REFERENCE  ELLIPSOID  USED 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 
c 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c  sfac  =  scale  factor  coorection 

c  hgd  =  horizontal  distance  c 
c  todis=  total  distance  in  traverse  c 
c  cegrid  =  grid  east  of  traverse  station  c 
c  cnqrid  =  grid  north  of  traverse  station  c 
c  difazi  =  azimuth  misclosure  c 
c  difdis  =  distance  misclosure  c 
c  coraz  =  angular  correction  per  station  c 
ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
C  c 
c  gridaz  =  subroutine  for  computing  azimuth  c 
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oooo  nnnnnnnnnnnnnn 


utm 


dmsr 


rdms 


cccccccco 


between  two  traverse  stations 
subroutine  for  computing  the  grid 
"Coordinates  from  known  distance 
and  azimuth 

subroutine  for  converting  the  angle 
from  degrees,  minutes,  and  seconds  to 
radians 

subroutine  for  converting  the  angle 
from  radians  to  degrees,  minutes,  and 
seconds 


ccccccccccccccccccccccccccccccccccccccco 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

C 

C 


C 

c 


c 

c 


DIMENSION  LAT(3) ,K(3) ,HGD(30) , HAGL( 30) , HANGf 30) 

DIMENSION  DHAGL(3o) ,mAAGL( 30) ’SHAGL(30) ,NAMt(3,30) , NAME1 ( 5 ) 
DIMENSION  NAME3( 2)  ,f^lAME4( 5)  ,MDlST( 30) 


* 

* 

* 

* 

* 

* 


VARIABLE  DECLARATION 


DOUBLE  PRECISION  AN90R,AN180R,AN360R.AZC .AZFIR.AZIMUT .AZLAS , 
CEGRld(  30  T,  CCiAZRC  30j ,C0$AZ ,CN6RID(30 j , 

distan,difAzi,deltax(30).d£ltay(30),dOmmyi, 

DUMMY2  DUMDIS  DIFDIS.TODiS 


DOUBLE  PRECISION 


DOUBLE  PRECISION 


* 

* 

* 

* 


DOUBLE  PRECISION 


NUM0 , NUM360; NUM1 . SUMDXJ  SUMDY^  STDD( 30 ) , STDA( 30) 
XGD,VGD,NUMf80,NL)M90,Wf*RED(3;? 

ANGS( 30) ,CFAZS(  30) ,DUM1D .DUM1M. DUM1S 
,  EAZIS ,  EAZ2S  .CORDxi ,  CORDYl , TEMS ,AZS21  ,AZS34 

SUMIX1 . SUMIX2 , SUMFY , SUNFX , SUMCOY . SUMCOX , 

ssta,s$td.sig){.sigy,sigxy.semaj  .Semin, 

SETA, SETAA.SSl6X.SSiGY.SSiGXY.DfGl, MINI. 

S ECOKl  1 ,  SEMAJ1 ,  SEMIN  1  ,  T!iMP01(  36)  , TEMP02(30) , 
DDIFFX  DDIFFY 


REAL*8  A.C,D,H,R,V,X,Y,GD.DH,ZD.CON.CUVaLAT,PHE,CUVRE.ELEV 

REAL*8  HOl^T ! Sdl5T!MC)I^Ti£SQft,GftAMMA1AN6^AN6VJ.Pi  , AGU ,  ^1 ,  F2 

REAL*8  K,KO.fl.E2.6l,Q2,ftR,VV  ECEN.FAC.SfAC.H6D 

REAL*8  F  CH6D.fiSDiST  DDM.HA.HAGL.CADISt.DG.fiMS.DS.DDANG.MMANG 

REAL*8  GNl.GNf .GN3,GN4,Gfl,6E2.Gf3.GE4.AZ2i,AZ34.6lST21,DIST34 

REAL*8  SI , $2 . D 1 , D2  S V , ANG  SHA6L , SANG , &D , MM , SS , HANG , RRR 

REAL*8  DUM1,&UM2,DUM3  DUM4  XXXX.YYYY 

REAL*8  IDD1  IMM1 , IDD2  IMM2,DVV,MVV, CUVRER.CUVR, CURVR 

REAL*8  ZINElKOSE 

INTEGER  I.N.IDl , ID2. IM1. IM2.DV,MV.DHAGL.MHAGL,DANG.MANG,UNIT 
INTEGER  ANG6(30).AN6M(3fl).A6Jl’ERft,AZD2i.AZM2i.AZD34,AZM34 
INTEGER  SSHAGLX36).CFAZDjr60LcfAZM(30J.MiN2,DE62, 

*  EAZ1D, EAZ lM , EAZ2  D , EA^2M , TEMD , ltMM 

DATA  DUM1/0.  30480061D0/ .DUM2/500000.  0D0/,DUM3/0. 000001D0/ 

DATA  DUM4/3600. 0D0/ , SIGANG/5. 0D0/.DUM5/0. 000005D0/,DUM6/0. 005D0/ 
DATA  TODIS/O.  0D0/,DUM10/0.  0000048481368D0/ 

DATA  NUMO/O. 0D0/.NUM360/360. 0D0/.NUM1/1.  0D0/.SUMDX/0.  0D0/, 

*  SUMDY/0. ODO/ , NUM180/180.  0D0/.NUM90/90.  OdO/ 


READ( 5,10 
READC5  15 
READC5  15 
READC5  16 
READ(5,16 
READC5  16 
READ(5  16 


AaF,K0 
l6l  IM1.S1 
ID2, IM2 ,S2 
GN1.GE1 
GN2  GE2 
GN3 ,GE3 
GN4  GE4 
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oooo  o  ooo  nnn  o  oo  o  oo  oo  o  oooooo 


READ( 5 , 1 1 )N , ELEV 

CALL  GRIDAZ  (GE1 , GN1 , GE2 , GN2 , AZ21 , DIST21 ) 
CALL  GRIDAZ  (GE3 |GN3 ,GE4 ,GN4 ,AZ34 ,DIST34) 

IDD1  =  FLOAT  C ID1 ) 

IMM1  =  FLOAT (IM1) 

CALL  DMSR  ( IDD1 . IMM1 ,S1 ,D1) 

LATfl)  =  6 1 
IDD2  =  FL0AT(ID2) 

I MM2  =  FLOAT U M2) 

CALL  DMSR  !IDD2.IMM2,S2,D2) 

LAT ( 2 )  =62 

LAT { 3 )  =  (Dl+D2)/2. 0D0 
ESQR  =  2.0D0*(1.0D07fJ-(1.0D0/F)**2 
X  =  A*DSQRT(1.  ODO-ESQR) 

Y  =  1.0D0-(ESQR*(DSIN(LAT(3)))**2) 

R  -  X/Y 


DETERMINATION  OF  THE  SCALE  FACTOR  FOR  UTM. 


El  =  DABS(DUM2-GE2) 

E2  =  DABS(DUM2-GE3) 

ECEN  =  ESQR/C1.  ODO-ESQR) 

01  =  DUM3*E1 
02  =  DUM3*E2 

QPRIME  =  ((Ql**2)+(Ql*Q2)+(Q2**2))/3.  0D0 
D01  M=1 ,3 

RR=A/( 1. 0D0-ESQR*(DSIN( LAT(M)))**2)**0.  5 

FI  =  ( 1. ODO+ECEN*DCOS( LAT(M) ) )*( 10.  0D0**12) 

F2  =  2. ODO*(RR**2)*(KO**2) 

FAC=  F1/F2 

C0NTKI®E=K0*(1'  0D0+FAC*QpRIME+(°-  00003D0*(QPRIME**2))) 
SFAC  =  K!3) 

SFAC  =  6.  0D0/((1.  0D0/K(l))+(4.  0D0/K(3) )+( 1.  0D0/K(2))) 


READ 

READ 


(5,20)  UNIT 

(5,18)(NAME1(L),L=1,5) 


DETERMINATION  OF  THE  HORIZONTAL  DISTANCES 
DO  1000  J=1 ,N 

READ( 5 , 14)(NAME( L, J) , L=1 , 3) 

READ! 5 , 12 )I ,DH,DV  MV  SV 
READf5.13)  $DI$T 
MDIST(J)  =  SDIST 

IN  CASE  OF  THE  LENGTH'S  UNIT  IS  IN  FEET,  THEN  CONVERSE  TO  METER 

IF(UNIT.  EQ.  1)  THEN 
SDIST  =  SDIST*DUM1 
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ooo  o  o  o  o  o  o  oooooo 


DH  =  DH*DUM1 
ELEV  =  ELEV*DUM1 
END  4F 


READ(  5^172DHAGL^J'>^MHAG|_^J I ,  SHAGL(J) 

MM  =  DFLOAT(MHAGL(J)) 

SS  =  SHAGL(J) 

CALL  DMSR  (  DD,MM,SS,RRR  ) 
HANGfJ)  =  RAR 
STDA(J)  =  SIGANG 


DETERMINATION  OF  HORIZONTAL  DISTANCES  WHEN  THE  DIFFERENT 
ELEVATION  IS  APPROXIMATELY  KNOWN. 

I F( I. EQ.  1)THEN 

IFfDH.NE.  O.OJTHEN 
DDM  =  SDIST 
IFQDDM.GE.  3300.0)THEN 
DO  100  KK=1 ,3 

CURV  =[0. 016192D0*DDM) 

CALL  DMSR  (NUM0,NUM0  ‘ 


IN 


_ ...  „  O.NUMO.CURV.CURVR) 

D=DSQRT(DDM**2-{DH*DC0S(CURVR)r*2)-DH*DSIN(CURVR) 


100 


CONTINUE 


DDM 


HDIST  = 

ELSE 

HDIST  =  DSQRT(DDM**2-DH**2) 

END  IF 

ELSE 

HDIST  =  SDIST 
END  IF 
END  IF 

DETERMINATION  OF  HORIZONTAL  DISTANCE  WHEN  ZENITH  DISTANCE  IS  KNOWN 
IF(I. EQ. 2)THEN 


DVV  =  FLOAT(DV; 


-  >  ') 

MW  =  FLOATf MV) 

CALL  DMSR  (DVV'mVV,SV,VANG) 

CALL  DMSR  (NUM§0 .NUMO . NUM0'AN90R) 
CALL  DMSR  ( NUM  18(5 ,  NUMd ,  NUM6  ,AN180R) 
ZD  =  AN90R-VANG 
DH  =  SDIST*DSIN(ZD) 

DS  =  SDIST*DCOSfZD) 

IFCSDIST.GE.  3300.01  THEN 
CUVRE  =  0.  0 1 39 25100* SD  1ST 
CUV  =  0. 016192D0*SDIST 
CALL  DMSR  (NUMO,NUMO,CUVRE,CUVRER) 
CALL  DMSR  (NUMO’NUMO,CUV,CUVR) 


C 

GRAMMA 
HDIST 
ELSE 
HDIST  = 
END  IF 
END  IF 


=  (AN90R-ZD)+CUVRER 
=  (AN90R+CUVR) 

=  AN180R-(C+PHE)  /  s 

=  ( SD I ST*DSIN( GRAMMA) )/DSIN(C) 

DS 
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oooo  o  o  o  o  o  o  oo  oo  ooo 


HH=  (ELEV+DH) 

determination^  the  horizontal  distance  on  geoid 

GD  =  R*HDIST/( R+H) 

DETERMINATION  OF  THE  GRID  DISTANCES 
HGD(J)  =  GD*SFAC 


TODIS  =  TODIS+HGDfJ) 

STDD(J)  =  HGD(J)*DUM5+DUM6 

WRITE(6,*)  'GRID  DISTANCE  =' ,HGD(J) 
WRITE(6!*;'  Y 


ELEV  =  HH 
1000  CONTINUE 

READ(5>19)(NAME3(L)>L=1,2) 
READ( 5 | 18) ( NAME4( L) ’ L=1 , 5 ) 


, DANG ,MANG, SANG 


CALL  OUTPUT! N, AZ21 .AZ34 .GN2 ,GE2 ,GN3 .GE3 ,HGD,DHAGL.MHAGL,SHAGL, 
*DANG. MANG ,  SANC5 ,  NAMEl .  NAME3 ,  AME4 ,  NAME  ,MdlST , GN1 ,  GEl , GN4 , GE4 ) 
dDANG  =  DFEOAT(D4NG) 

MMANG  =  D FLOAT (MANG) 

CALL  DMSR  (DDANG, MMANG, SANG, RRR  ) 

NN  =  N+l 
HANG(NN)  =  RRR 
STDA(NN)  =  SIGANG 


C 


C 


C 

C 

C 

C 

C 

c 


CALL 

CALL 


DMSR 

DMSR 


NUM360 ,NUM0 ,NUM0 ,AN360R  ) 
NUM180 , NUMO ,NUM0 ,AN180R  ) 


DO  2000  ADJ1  =  1,2 
XGD  =  GE2 
YGD  =  GN2 
AZFIR  =  AZ21 
DO  200  I  =  1  ,N 

AZIMUT  =  AZFIR+HANG( I) 

IF  (AZIMUT.  GE.AN360R)  THEN 
AZIMUT  =  AZIMUT-AN360R 
END  IF 

IF(AZIMUT.  GE.  AN180R)  THEN 
AZIMUT  =  AZIMUT-AN180R 
ELSE 

AZIMUT  =  AZIMUT+AN180R 
END  IF 

COAZR(I)  =  AZIMUT 


CALL  RDMS  (  AZIMUT,DUM1D,DUM1M,DUM1S  ) 
CFAZD(I)  =  idlNT(DUMlD) 

CFAZM(I)  =  IDINT(DUMIM) 

CFAZS(I)  =  DUM1S 


DISTAN  =  HGD(I) 

CALL  UTM  (  XGD, YGD, DISTAN, AZIMUT, DUMMY1 ,DUMMY2  ) 
WRITE(6,*)'I  =' ,1 
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ooo  ooo  o  o  oo  o  o  o  oooooooo 


WRITEC6 ,*) ' GRID  E=' . DUMMY 1. 
DELTAX(I)  =  DUMMYl-XGd 


GRID  N=‘ ,DUMMY2 


_  DELTAYCI)  =  DUMMY2-YGD 

17RITE(6  *j  DEPARTure  ='  DELTAX(I) , 1  LATITUde  =' 
ZlNE  =  DSINC AZIMUt) 

KOSE  =  DCOS(AZIMUT) 

SUMDX  =  SUMDX+DELTAX( I ) 

SUMDY  =  SUMDY+DELTAY( I ) 

CEGRID(I)  =  DUMMY1 
CNGRID(I)  =  DUMMY2 

CALL  GRIDAZ  (  XGD , YGD , DUMMY 1 , DUMMY2 , AZC , DUMDIS  ) 

CALL  RDMS  (  AZC.DUMID^UMIM^DUMIS  ) 

CBAZD(I)  =  IDINT(DUMID) 

CBAZM(I)  =  IDINT(DUMIM) 

CBAZS(I)  =  DUM1S 


IF  (ADJ1.  EQ.  1)  THEN 

TEMPOl(I)  =  CEGRID(I) 
TEMP02(I;  =  CNGRIDC I ) 
END  IF 


AZFIR  =  AZC 
XGD  =  DUMMY 1 
YGD  =  DUMMY2 
200  CONTINUE 

IF  (ADJ1.NE.  1)  THEN 

CORDX1  =  (  SUMDX 
CORDY1  =  l  SUMDY 
XGD  =  GE2 
YGD  =  GN2 


=  (  SUMDX-(GE3-GE2)  )/TODIS 
=  ^SUMDY-(GN3-GN2)  )/T0DIS 


DO  230  1=1, N 

DISTAN  =  HGD(I) 

XXXX  =  C0RDX1*HGD( I) 

YYYY  =  CORDYl*HGDa) 

DELTAX(I)  =  DEllTAX(I) 
DELTAY(I)  =  DELTAY(I) 


XXXX 

YYYY 


CEGRID(I)  = 
CNGRIDC I )  = 
DDIFFX  = 
DDIFFY 


XGD+  DELTAXfl' 

YGD+  DELTAYfl' 
CEGRID(I)-TEMP01(I) 
CNGRIDC I )-TEMP02( I ) 


C 

c 


230 


XGD 

YGD 

CONTINUE 


=  CEGRIDf I' 
=  CNGRIDC I! 


END  IF 

AZLAS  =  AZFIR+HANG(NN) 

IF  (  AZLAS.GE. AN360R  )  THEN 
AZLAS  =  AZLAS-AN360R 
END  IF 

IF  (AZLAS.  GT.AN180R)  THEN 


,DELTAY( I) 
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ooooo 


c 

c 


AZLAS 
ELSE 
AZLAS 
END  IF 


AZLAS-AN180R 

AZLAS+AN180R 


CALL  RDMS  (  AZLAS ,DUM1D.DUM1M,DUM1S  ) 
TEMD  =  IDINt(DUMlD) 

1 


TEMM  =  IDINT(DUM1M 
TEMS  =  DUM1S 


DIFAZI  =  AZLAS-AZ34 
WPRED(l)  =  DIFAZI  ,  N 

CORAZ  =  DI FAZ I/DFLOAT ( NN ) 

CALL  RDMS  (  DIFAZI ,DUM1D,DUM1M,DUM1S  ) 
EAZ1D  =  IDINt(DABS(  DUMlDn 
EAZ1M  =  IDINTC DABS(DUMIM)) 

EAZ1S  =  DABS(DUMIS) 

CALL  RDMS  (  CORAZ .DUM1D.DUM1M.DUM1S  ) 
EAZ2D  =  IDINT(DAB$(DUM1D)) 

EAZ2M  =  IDINTCDABS(DUMIM)) 

EAZ2S  =  DABS(DUMIS) 


IFfADJl.  EC 
WPREDC ‘ 
WPRED 
DIFDI. 
ERR 
END  IF 


1)THEN 

=  SUMDX  -  (  GE3-GE2  ) 

=  SUMDY  -  I  GN3-GN2  ) 

=  DSQRTf  WPRED(2)**2  +  WPRED(3)**2 
=  IDINT(  TODIS  /  DIFDIS  ) 


PRINT  RESULTS  OF  TRAVERSE  COMPUTATION 


CALL  RDMS  (AZ21 ,DUM1D.DUM1M,DUM1S) 
AZD21  =  iDINTtfiUMID) 

AZM21  =  IDINT(DUMIM) 

AZS21  =  DUM1S 

CALL  RDMS  (AZ34,DUM1DaDUM1M,DUM1S) 
AZD34  =  iDINT(d)UMlD) 

AZM34  =  IDINT(DUMIM) 

AZS34  =  DUM1S 


IF  (ADJ1.  EQ.  1) 
WRITEC8 ,29) 
WRITE(8  27^ 
WRITE*  8*30 
WRITE(8  31 
ELSE 

WRITE(8,28 
WRITE(8  27 
WRITE*  8,30 
WRITE(8  40 
END  IF 
WRITEC8 , 32 
WRITE* 8 , 33 
WRITE(8  30 
WRITE*  8,34 
WRITE*  8 ’ 35 


THEN 


AZD21 ,AZM21 ,AZS21 
GE2 ,  G^I2 


DO  240  1=1, N 
II  =  1+1 

CALL  RDMS(HANG( I) ,DUM1D,DUM1M.DUM1S) 
ANGD(I)  =  IDINT(DUMID) 

ANGM(I)  =  IDINT(DUMIM) 

ANGS(I)  =  DUM1S 

SSHAGL(i)  =  IDI NT ( SHAGL( I ) ) 
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IF  (ADJ1.  EQ.  1)  THEN 

. w 


WRITE( 8 , 


)DHAGL 
CFAZM( I 


'WflGD(I) 


C 

c 


240 


ELSE 

WRITE(8 , 37)ANGD 
CFAZS( 

END  IF 

WRITE(8 , 38)  II ,CEGRID( I) ,CNGRID( I) 
CONTINUE 

CALL  RDMS  (  HANG(NN),DUM1D,DUM1M,DUM1S) 
IF  (ADJ1.EQ.  1)  THEfl  ' 


®tWs^?fL(I)'CFAZD<I)' 

, ANGM( I ) , ANGS (I) ,CFAZD(I) ,C FAZM( I ) , 


DHAGL( 

MHAGLf 


=  DANG 
=  MANG 


ELSE 


SSHAGL(NN)  =  IDINT(SANG) 

_WRITE(8,36)DHAGL(NN),MHAGL(NN),SSHAGL(NN)  ,TEMD,TEMM  JEMS 
CALL  RDMS  (HANG(NN) ,DUM1D,DUM1M,DUM1S) 


ANGD(NN)  =  IOINTIDUMID, 

ANGM(NN)  =  IDINT(DUMIM) 

ANGS(NN)  =  DUM1S 

end  WR.ITE(8,37)ANGD(NN),ANGM(NN),ANGS(NN),TEMD  ,TEMM  JEMS 
WRITE(8,30) 

WRITEf8,39)AZp34,AZM34,AZS34,GE3,GN3 
WRITEC8  '  r 


*  \  i 

iT(A5ji.  EQ.  1)  THEN 

WRITE(8 , 4I)EAZ1D , EAZ1M , EAZ1S 
WRITE! 8 , 42 ) EAZ2D , EAZ2M , EAZ2S 
WRITE(8  43)TODIS 
WRITE!8,44)DIFDIS 
WRITE(8,45)  ERR 
END  IF 

CORRECTED  OBSERVED  ANGLES 

DO  250  1=1, NN 

HANG(I)  =  HANG(I)-CORAZ 
250  CONTINUE 

SUMDX  =  0. 0D0 
SUMDY  =  0. 0D0 


C 

c, 

c‘ 

c 


2000  CONTINUE 


SUMDX  =  GE3-CEGRID(N) 

SUMDY  =  GN3-CNGRIDCN) 

DIFDIS=  DSQRT  (  SUMDX**2  +  SUMDY**2  ) 

WRITEf 8,41)  EAZ1D . EAZ1M,EAZ1S 
WRITE(8’44)  DIFDlS 

10  FORMAT! 5X.F15.  5 . 5X , F15, 10 , 15X , F12.  8) 

11  FORMAT{26* ,I4,3$X.f7. 3)  ' 

12  FORMAT! 20 X ’ Il’4X,f9.3,llX,2I3,F10.4) 

13  FORMAT! 25X! Fl$.  6) 

14  FORMAT! 6X ,3A4) 

15  FORMAT(5X!2I3J14.8) 

16  FORMAT! 10*, Fl$.  5.15X.F15.  5) 

17  FORMATf 25X*2I3, FlO.  4) 

18  FORMAT! 5A4) 

19  FORMAT! 2A4.17X, 213, F10. 4) 

20  FORMAT;  II) 

' 27  FORMAT( 16X , 1 ::::::::::::::::::::::::::::: 1 ,///) 
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28 

29 

30 

31 
40 

32 

33 

34 

35 

36 

37 

38 
C39 
C41 
C 

C42 

C 

C43 

C44 

C45 


FORMAT(16X, 'ADJUSTED  TRAVERSE  COMPUTATION1 ) 
FORMAT? 15X  'UNADJUSTED  TRAVERSE  COMPUTATION') 

cnDMATr  i - - - - - 


FORMAT? 

* _ _  I  N 

FORMAT? 1 j  5TN 
MAT?-'  j  ^TN 


rur\r 
*. ) 
FOfo 
*TEi 
FOR 


flAT 


( 


FORMAT? 1 ; 
*N0RTH  : 1 
'I 


FORMAT? 
FORMAT? 
FORMAT? 1 

*  Jrr 


o. 


i  i 
1 1  ,3X 


OBSERVED  :  FORWARD 

CORR.  OBS. :  FORWARD 

HOR.  ANG.  :  AZIMUTHS 

D  M  S  :  D  M  S 

'  ,13X,I3.1X,I2,1X.F5.2,34X.  * 
'  36X’Fld.3’lX.Fll.3.1*  'I') 

, r: Y,i3,lX,l2,iX,I2,4x,i3,lX, 


GRID 
GRID 
DIST. 

(M.  ) 

'I') 

12 , IX , F5.  2 , IX , F8. 3,25X 


COORDINATE  (M 
ADJUSTED  COORDINA 

GRID  EAST  :  GRID 


FbRMATf 'I' ,3X,':  ' ,13, IX, 12, IX, F5.  2, IX, 13, IX, 12, IX, F5.  2,1X,F8.  3, 
*25X  1 1 1 ) 

FORMAT? ' 1 '  1 3 , ' : '  36X, F10. 3 . IX . FI 1. 3 , IX  ' I' ) 

FORMAT(lX/KNdWN  tiATA' ,  7X  ,  li ,  1>? ,  12.1)? .  F$.  2 , 11X  ,  F10.  3 ,  IX  ,  Fll.  3.//) 
FORMAjr ANGULAR  ERROR  ='  ,2X ,  13 , IX  ' D'  IX ,  12  IX  ' M',1^ , F5. 


»?  1  v 

tormAt? 

*2, IX. 'S 

f6rmAt? 

FORMAT? 

FORMAT? 

STOP 

END 


NGULAR  CORR./  STATION  =' ,2X , 13 , IX , ' D',1X , 12 , IX , ' M'.IX , F5. 


loTAL  DISTANCES 
LINEAR  ERROR 
LINEAR  CLOSURE 


=' ,2X,F10. 3) 

='  2X  F10. 3) 

='  3X  ' 1  r  16,///////////) 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

C  SUBROUTINE  TO  CALCULATE  GRID  AZIMUTH  AND  DISTANCE  BETWEEN  C 
C  TWO  KNOWN  STATIONS  WITH  THEIRS  GRID  NORTH  AND  GRID  EAST.  C 

c  c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

SUBROUTINE  GRIDAZ  (Xl.Yl.XZ.YZ.ANGLR.DISTIZ) 

c 

c 

REAL*8  XI ,X2 . Y1 , Y2 .ANGLR , ANGO ,ANG90 ,ANG180 .ANG270 ,DIFX ,DIFY .DIST12 
REAL*8  AN(j90ft,AN18(!)R ,AN2^0R 


C 

c 

c 

c 


DATA  ANGO/O.  ODO/ .ANG90/90.  0D0/.ANG180/180.  ODO/ .ANG270/270.  ODO/ 


CALL  DMSR 
CALL  DMSR 
CALL  DMSR 

DIFX  =  X2-X1 
DIFY  =  Y2-Y1 


ANG90 , ANGO .ANGO .ANG90R  ) 
ANG18?)  ,ANG()  ,ANG()  ,AN180R 
ANG270, ANGO  ANGO, AN270R 


C 

C 


IF? (DIFX.  EQ.  ANGO).  AND.  (DIFY.  EQ.  ANGO))  THEN 
ANGLR  =  ANGO 

ELSE  IFfDIFX.  EQ.ANGO)  THEN 
IF(DIFY.GT.ANGO)  THEN 
ANGLR  =  ANGO 
END  IF 

ANGLR  =  AN180R 

ELSE 

IF(DIFX.GT.  ANGO)  THEN 

ELS^GLR  =  ANG90R-DATAN(DIFY/DIFX) 

ANGLR  =  AN270R-DATAN(DIFY/DIFX) 
END  IF 

END  IF 

DIST12  =  DSQRT(DIFX**2+DIFY**2) 
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RETURN 

END 


0000000000000000000000000000000000000000000000000000000000 
0  (3 

0  SUBROUTINE  OUTPUT  TO  PRINT  OUT  THE  RESULT  FOR  USING  IN  0 

0  THE  NEXT  COMPUTATION.  0 

0  0 

(3  (3  (3  (3(3  (3(3  (3  (3  (3(3  (3  (3  (3(3(3  (3  (3(3  (3(30(3  (3(90(3(3  (3(3  (3(3(3(3  (30(3(3(3  (3@0(3(3(3@P(3(9(3@(3(3(3(300@ 


c 

c 


c 

c 


10 


c 

c 


SUBROUTINE  OUTPUT ( N . A , B , G2 , E2 . G3 . E3 , RD , DH , MH , SH , DA , MA , SA , NAME1 , 
*NAME3 , NAME4, NAME, MD f Sf , dl ,  £l  ,d4 ,  E:4) 


G,G2.G 


DH(N),MH(N),SH(N) ,NAME(3,N) , NAME 1 ( 5 ) , NAME3 ( 2 ) 


DIMENSION  RD 
\NAME4(5).MD:  .  „ 

kEAL*8V  A\.&,DG,G2aG3.E2.E3  .RD.AA.BB.HH.SH ,SA,SSA,SIGA, SIGD 
REAL*8  M6l§T ,MDI§ , RAD , 64 , £l , ^4 , CUMl 
INTEGER  N,DH,MH,DA,MA 

DATA  SIGA/5. 0D0/,DUMl/0. 30480061D0/ 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 


=  fl+i 


,240' 

250 

260 

250 

270 

280 

290 

280 

240 

,280 

300 

,310 

310 

300 

280 

240 


'NAME 

'NAME 

'NAME 

NAME4 


WR 

WR 

WR 

WR 

WR 

WR 

WR 

WR 

WR 

WR 

WR 


TE( 

TE( 

TE( 

TE< 

TE 

TE 

TE 

TE? 

TE( 

TE 

TE 


,230' 
,100 
,110 
120 
110 
125 
,140 
130 
140 
,100 
,140 


N1 

=  N+l 

DO 

10 

C\J 

II 

l — l 

WRITE! 6 , 

* 

CONI 

r1-- 

WR] 

[TE( 

'6,150 

WR] 

[TE( 

6  150 

WR] 

[TE( 

6  140 

WR] 

[TE( 

6  100 

\y 

)(NAME( 


WRITE(6,230) 
WRITE(9 1 170) 
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c 

20 

C 

C 

40 

100 

110 

120 

125 

130 

140 

150 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 


WRI 

WRI 

WRI 

WRI 

WRI 

WRI 

WRI 

WRI 

WRI 


'9,180' 
'9,490 
'9,180 
'9  170 
'9  200 
'9  210 
'9  200 
'9  170 


E(9,200; 


DO  20  J=1,N-1 

MDIS  =  MDIST! J)*DUM1 

RRD  =  RD! J)/DUM1 

SIGD  =  (RD!J)*0. 00001D0)+0. 01D0 

*  ;2R2D? iSOT . 2 ’ , ( NAME( L ,  ) ,  L= ! , 2 ) , MDI S , 

CONTINUE 

MDIS  =  MDIST(N)*DUM1 

RRD  =  RD!N)/DUM1  x  % 

ytWRITE(9^220)(NAME(L,N),L=l,2),(NAME3(L),L=l,2),MDIS,MDIST(N) 

$IGD  =  (RD!N)*0. 00001D0)+0.  01D0 
WRITE(9,200) 

WRITE(9 , 170) 


FORMAT(2I3 . F6.  2) 

F0RMATT9X, f****^**** 

******>T \ 

FORMAT^X,'* 

FORMAT^,'*  MEASURED  HORIZONTAL  ANGLE  * 


********************************************** 

*  *  * 

*  * 


FORMAT! 9^ , 1  *******************************  DEGREES  *  MINUTES  *  SEC 
*ONDS  *r) 

FORMAT^X,'*  AT  *  TO 

FORMAT^,'*  *  *  *  * 

FORMAT !  9^ . 1  * 1 , 3X . 2A4 , 3X , 1  * 1 ,3X,2A4,3X, ,3X,I3,3X, ,4X,I2,3X, 
2X,F$,2,2X 

formAtc  Sx . 1  kttkmhmmmmMmttmmmmmmmmm 

FORMAT! 5X,‘#  »  u 

*  ,  .  #') 

FORMAT! 5X'#  DISTANCES 

*D  DISTANCES  #r ) 

FORMAT!  5X/#  ,  #  # 

*  #  #  ) 

FORMAT! 5X  '#  FROM  #  TO  # 

*  .# - Wl...  J\) 


# 

# 


MEASURED  DISTANCES 
# 

METERS  #  FEET 


# 

» 

# 


REDUCE 

METERS 


FORMAT! 5X .2A4: 2A4, '# ' ,2X , F8.  3 , 2X , '# 1 ,2X , F8. 2 , 2X , 1  § '  , 2X , 
*F8. 3 ,  IX ,  #  , £x , F&. 2 ,CX ,  A  ) 

FORMAt('i') 

F0RMAT(9X  ^  ******************************************************* 

FORMAT! 9X , 1  *  *  *' 

FORMAT! 9X,1*  NAME  *  COORDINATES  *' 

FORMAT! 9X , 1  *  OF  *********************************. 

FORMAT! 9X, 1  *  *  *  *' 

FORMAT ( 9X , 1  *  KNOWN  POSITIONS  *  GRID  NORTH  *  GRID  EAST  *' 

F^RMAT( 9X , ' * 1 ,5A4  . 2X . FI  1.  3 ,2X, ' * '  2X . Fll. 3 , 2X, ' * 1 ) 

FORMAT! 9X  '  *'  4X , 3A4 , 4x ,  ut , 2X  Fll.  3 , Cx , ' * '  ,2X  Fll. 3 ,2X , 1  * ' ) 

RETURN 
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END 


(3(9(9(3(3(3(3(3(3(3Q(3(3PP(3(3(3(3(3(3(3(3(3(B(9(3(3PP(3(BPP 
0  @ 

0  CONVERT  DEGREES,  MINUTES,  0 

0  AND  SECONDS  TO  &ADIANS  0 

0  0 

0000000000000000000000000000000000 


SUBROUTINE  DMSR  (  D1,M1,S1,R1  ) 

REAL*8  D1,M1,S1,R1,PI1 ,N60 ,N3600 ,N180 

DATA  N60/60. 0D0/.N3600/3600. 0D0/,N180/180.  0D0/ 

PI1  =  4.  0D0*DATAN( 1. 0D0) 

R1  =  ((D1+(M1/N60)+(S17N3600))*PI1)/N180 

RETURN 

END 


C 

C 

c 


c 

c 


000000000000000000000000000000000 
@  @ 

0  COMPUTE  THE  COORDINATES  0 

0  FROM  ONE  KNOWN  POINT  WITH  0 

0  MEASURED  DISTANCE  AND  ANGLE  0 

0  0 

000000000000000000000000000000000 

SUBROUTINE  UTM  (XI , Y1 , DIS12 .ANG12R, X2 , Y2  ) 

REAL*8  XI ,Y1 ,X2,Y2 ,DIS12,ANG12R,DIFX ,DIFY 

DIFX  =  DIS12*(DSIN(ANG12R) ) 

DIFY  =  DIS12*(DCOS(ANG12R)) 

X2  =  Xl+DIFX 
Y2  =  Yl+DIFY 

RETURN 

END 

000000000000000000000000000000000 
0  0 

@  CONVERT  ANGLE  FROM  RADIAN  0 

0  TO  DEGREE, MINUTES, SECONDS  0 

0  0 

000000000000000000000000000000000 

SUBROUTINE  RDMS  (  RR1,D2,M2,S2  ) 

REAL*8  RR1 ,D2 ,M2 ,S2 , PI1 , N60 , N180 ,TDEG,DIF,TMIN 

DATA  N60/60. 0D0/,N180/180. 0D0/ 

PI1  =  4. 0D0*DATAN(1. 0D0) 

TDEG=  f N180*RR1)/PI1 
D2  =  DFLOAT( IDINT(TDEG)) 

DIF  =  TDEG-D2 
TMIN=  DI F*N60 
M2  =  DFLOATflDINT(TMIN)) 

DIF  =  TMIN-M2 
S2  =  DIF*N60 

RETURN 

END 
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APPENDIX  C 

INDTRA  FORTRAN  PROGRAM 


INDTRA  Fortran  program  is  used  for  computing  and  adjusting  traverse  station 
position  by  least  squares  method  of  observation  equations  (indirect  observation)  and 
simultaneously  adjust  the  intersection  points  which  are  observed  at  each  traverse 
station. 


(3(3  (3(3  (3  (3  (3  (3  (3  (3(3  (3  @  (3  (3  @  (B  (3(3@(3  (3(3@(3  (3  (3  (3  (30(3 
00  (30 

00  LEAST  SQUARE  ADJUSTMENT  OF  00 

00  OBSERVATION  EQUATIONS  METHOD  (3(3 

(3(3  (3(3 

000000000000000000000000000000000000000 


000000000000000000000000000000000000000000000000000000 


(3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

N  ,NA  =  NUMBER  OF  OBS.  ANGLES  0 

Nftl  =  NUMBER  OF  RESECTION  PT.  0 

ND  =  OBS.  DISTANCES  0 

NT  =  OBS.  DIST  +  OBS.  ANG.  0 

NC  =  TRAV.  STN.  +  RESECT.  PT.  0 

(NOT  INCLUDE  KNOWN  STN. )  0 

NZ1  =  NUMBER  OF  RESECT.  ANGLE  AT  EACH  TRAV.  STN.  0 
EGRID  =  GRID  EAST  OF  KNOWN  POINT  0 

NGRID  =  GRID  NORTH  OF  KNOWN  POINT  0 

D.MaS  =  READ  ANGLE  IN  DEGREES,  MINUTES,  SECONDS  0 
RAN6.  DIST  =  GRID  DISTANCES  0 

NCO, £C0  =  GRID  NORTH  AND  EAST  OF  INTERSECTION  PT.  0 
STDA  =  STANDARD  DEVIATION  OF  ANG.  0 

STDD  =  STANDARD  DEVIATION  OF  DIST.  0 

NP2  =  NUMBER  OF  UNKNOWN  DX , DY  0 

CEGRID.CNGRID  =  TRAV.  STN.  CfiORDI.  0 

INDLSQ  =  SUBROUTINE  ADJUSTING  THE  COORDINATES  0 

BY  OBSERVATION  EQUATIONS  0 

CALFM  =  SUBROUTINE  FOR  COMPUTING  CONSTANT  VECTOR  0 
CALAM  =  SUBROUTINE  FOR  COMPUTING  THE  COEFFICIENTS  0 


OF  UNKNOWN  PARAMETERS 


0 

0 


(3(3  (3  (3  (3  (3(3  (3(3  (3(3  (3  (3  (3  (3  (3(3  (3(3  (3  (3  (3  (3  (3  (3  (3  (3  (3  (3(3  (3(3  (3  (3  (3  (3  (3  (3  (3  (3  (3  (3  (3  (3(3  (3(3  (3(3  (3  (3  (3(3(3 


* 

* 

* 

* 

* 

* 

I 

* 

* 

* 

[ 

* 

* 


DOUBLE  PRECISION 


DOUBLE  PRECISION 


DOUBLE  PRECISION 


NGRID(4).EGRID(4J.AZ21 .DIST21.AZ34, 
DIST34,X(jD,  YGD  !/az  f- 1  r,nUm360,  nUmo  , 

AN360R ! DD .MM , S$ , R , AN6l(40) . $TDA(40) . 

dumio , &AN6,DiST(4d) ,ncu,ec6.cngrid(4o) , 

CEGRlflUO) ,STDD(40) ,AZIMUT,DlS(40) . 
coazr(4o) j Dummy 1.DUMMY2 ,azC,cfazs(4o) , 
DUM1D ,DUMiM,DUMl$ vDELX(40) ,DELY(4u1 
SUMDX ,  SUMDY,  CBAZSj  40 ) ,  TTCOk ( 40 j  , TTCOYjf  40 ) , 
AZLAS.DIFAZf ,WPRED(3) . EAZ1S XORAZ . DIFDIS , 
EAZ2S ,TODIS ,ASNY(40) ,ASEX(4D) , DUMi 1 , 

DUM12 

WMf  80 1 , 80 ). STD0(80) , FMf 80 ) . BMC  80 '  80) . 

BMT(8d ,80 j ,BMTWM(od.8U) ,NM(80A8d) ,TM(80) , 
DELTA(DO) ,NMI(80,80) ,WKlO(428D) ,VM( 1 ,80) , 
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c 

c 


c 

c 


c 


* 

ye 

* 


INTEGER 


STDAD,XC0X(80) ,YCOY(80) 

NT.NA,NP,NP2.ND,NI .PN.ASP, PRINTD.N.NSTA, 

D .M.S . CF/ZDf  40)  ,CF/4ZM(40) ’CBAZDf  40  j ,  K ,  Kl(, 
k4,£bAZM(40),EA21D.EAZ1M,£AZ2D,EAZ2M2, 

NUM(40)  ,ERR,NRl,NZi(40) ,  K1 0  ,  KOKl ,  1 1  ,lil,KlC 


DATA  NUM360/360.  0D0/,NUM0/0. 0D0/ .DUM10/2. 0D0/ 
DATA  DUM11/0.  00001D0ADUM12/0.  0260/ , SUMDX/0.  0D0/ 
DATA  SUMDY/0.  ODO/,TODiS/O.ODO/ 

DATA  PW/0/ , PRINTD/1/ 

PW  =  1  UNWEIGHT 

PRINTD  =  0  NOT  PRINT  DETAILS  OF  MATRIX 


DETERMINE  BACK  AZIMUT  AT  STARTING  STATION  AND  FORWARD 
AZIMUTH  AT  CLOSING  STATION 


DO  1  1=1,4 

READ (5 , *) 

i  NGRIDl 

READ(5  * 

i  EGRIDI 

CONTINUE 

CALL  GRIDAZ  i 

f EGRIDC 

CALL  GRIDAZ  i 

CEGRIDC 

(1) 

!) ,NGRID(2) ,EGRID(1) ,NGRID( 1) , 

I)!NGRID(3)!EGRID(4)!NGRID(4)! 


AZ21.DIST21) 

AZ34.DIST34) 


COMPUTE  APPROXIMATED  COORDINATE  OF  TRAVERSE  STATION  POSITIONS 


coot? 

c 

c 

coooo 

c 

c 

c 


c 

c 

c 

c 

c 

c 


READ  NUMBER  OF  OBSERVED  ANGLES 
READ(5,*)  N 

READ  NUMBER  OF  RESECTION  POINT 
READ(5 ,*)  NR1 

NSTA  =  N-l 
XGD  =  EGRID(2) 

YGD  =  NGRID(2) 

AZFIR  =  AZ21 

CALL  DMSR  (NUM360.NUM0 ,NUM0 .AN360R) 

K  =  0 

DO  100  1=1, NSTA 

READ  ANGt!E  IN  DEGREE,  MINUTE,  SECOND 

read(dV^f&)aVo; 

MM  =  DFLOATfM 
SS  =  DFLOATfS' 


CALL  DMSR  (DD,MM,SS,R) 

ANGl(I)  =  R 
STDA(I)  =  DUM10 

READ  REDUCE  GRID  DISTANCES 
READ(5  *)  RANG 

DfST( I )  =  RANG 
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c 

c 

c 

c 

c 


READ(5,*)  NZ1(I) 


IF  (RANG.  EQ.  0.  0)  THEN 

READ  APPROX.  GRID  COORDINATES  OF  RESECTION  STATION 
READ( 5 ,*)  NCO 
READ(5  *)  ECO 


END  IF 


CNGRID(I)  =  NCO 
CEGRID(I)  =  ECO 
GO  TO  100 


K  =  K+l 

DIS(K)  =  DIST(I) 

STDD(K)  =  (DIS(K)*DUM11)  +  DUM12 

AZIMUT  =  AZFIR  +  ANGl(I) 

IF  (AZIMUT. GT.AN360R)  THEN 
AZIMUT  =  AZIMUT  -  AN360R 
END  IF 


C 

cooooo 

c 


cooooo 

c 


COAZR(K)  =  AZIMUT 
FOR  PRINTING 


CALL  RDMS  (AZIMUT,DUM1DaDUM1M,DUM1S) 
CFAZD(K)  =  1 DINT(OUMID) 
CFAZM(K)  =  IDINT(DUMIM) 
CFAZS(K)  =  DUM1S 


CALL  UTM  (XGD.YGD, RANG, AZIMUT, DUMMY!, DUMMY2) 


C 

C 

C 


c 


DELX(K)  = 
DELY(K)  = 


SUMDX 

SUMDY 

CEGRID 

CNGRID 


DUMMY1-XGD 
DUMMY2-YGD 

=  SUMDX  +  DELX(K) 
=  SUMDY  +  DELY(K) 


(1): 


DUMMY1 

DUMMY2 


C 

C 


CALL  GRIDAZ  (DUMMY1.DUMMY2, XGD.YGD, AZC.DUM1S) 
FOR  PRINTING 

CALL  RDMS  (AZC,DUM1D,DUM1M,DUM1S) 

CBAZD(K)  =  IDINT(DUMID) 

CBAZM(K)  =  IDINT(DUMIM) 

CBAZS(K)  =  DUM1S 


AZFIR  =  AZC 
XCOX(K)  =  DUMMY 1 
YCOY(K)  =  DUMMY2 

XGD  =  DUMMY1 

YGD  =  DUMMY2 

TODIS  =  TODIS  +  RANGE 


C 

C 

C 

C 

100  CONTINUE 
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ooo 


c 

c 

c 


K10  =  K^l 


9 

C 

C 

C 


c 


c 

c 

c 


c 


c 


c 

C@(3@(9@ 

C 

c 


c 

c 


C@(3(3(3(3 

c 


c 


101 

c 

c 


DO  9  I  =1 , K10 

rrcoxm  =  xcoxrr 
TTCOY(I)  =  YC0Y(I' 
CONTINUE 


KK  =  K+l 

READ  (5.105)  DAMAS 
D  t)  =  DFLdAtfD) 

MM  =  DFLOAT  M) 

SS  =  DFLOAT(S) 

READ  (5  *)  RANG 

DiST(N)  =  RANG 

CALL  DMSR  (DD,MM,SS,R) 

ANG1(N)  =  R 
STDA(N)  =  DUM10 

AZLAS  =  AZFIR  +  ANG1(N) 

IF  (AZLAS.GT.  AN360R)  THEN 
AZLAS  =  AZLAS  -  AN360R 
END  IF 

DIFAZI  =  AZLAS  -  AZ34 

WPRED(1)=  DIFAZI 

CORAZ  =  DIFAZI  /  DFLOAT(KK) 

FOR  PRINTING 

CALL  RDMS  (DIFAZI ,DUM1D,DUM1M,DUM1S) 

EAZ1D  =  I D I NT(  ABS  (DUM1D)  ) 
EAZ1M  =  IDINTC  ABS  (DUM1M)  ) 
EAZ1S  =  ABS  (DUM1S) 

CALL  RDMS  (C0RAZ,DUM1D,DUM1M,DUM1S) 

EAZ2D  =  I D I NT(  ABS  (DUM1D)  ) 
EAZ2M  =  IDINTf  ABS  (DUM1M)  ) 
EAZ2S  =  ABS  (DUM1S) 


WPRED(2j  = 


WPRED(3; 

DIFDIS 

ERR 


SUMDX  -  (EGRID(3)  -  EGRID(2) ) 
SUMDY  -  (NGRID(3)  “  NGRID(2) ) 


=  DSQRT (  (WPRED(2))**2  +  (WPRED(3))**2  ) 
=  I  DINT  L  TODIS  >  DIFDIS  ) 


READ  NUMBER  OF  MEASURED  ANGLES  AT  EACH  TRAVERSE  STATION 


DO  101  J=1 , KK 
READ (5,*)  N2 
nl)m(J)  = 

CONTINUE 


N2 


DO  108  I  =  1,NR1 
II  =  K10+I 
READ(5,*)  DUMMY1 
READ(5!*j  DUMMY2 
TTCOY(ll)  =  DUMMY1 
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108 

C 

C 

C 

C 

m 

m 

c 


c 

102 

C 

m 

c 

c 

C@@ 

c 

c 

c 

c 

C@@ 

c 

c 

C@@ 

C 

C 

C 

C@(3 

C 

C@(3 

C 

C 


TTCOX(II)  =  DUMMY2 
CONTINUE  ~ 

ASP  =  NSTA  -  1 

KEPP  APPROX.  GRID  COORDINATES  OF  EACH  STATION  IN 
ASNY  AND  ASEX 


DO  102  1=1, ASP 

ASNY(I)  =  CNGRID(I) 

ASEX(I)  =  CEGRID(I) 

CONTINUE 

NUMBER  OF  OBSERVED  ANGLES  PLUS  OBSERVED  DISTANCES 
NT  =  N+K 

NUMBER  OF  OBSERVED  ANGLES 
NA  =  N 

NUMBER  OF  OBSERVED  DISTANCES 
ND  =  K 

NUMBER  OF  TRAVERSE  STATIONS  NOT  INCLUDED  KNOWN  STATIONS 
NP  =  ASP 

NUMBER  OF  UNKNOWN  DX  AND  DY 
NP2  =  K10*2  +  NR1*2 

NI  =  NP2**2  +  3*NP2 

NUMBER  OF  TRAVERSE  STATION  (NOT  INCLUDE  KNOWN  STATION) 
PLUSE  NUMBER  OF  INTERSECTION  POINT 

NC  =  K10+NR1 


C(3@(3(9(3  CALL  SUBROUTINE  TO  ADJUST  STATION  POSITIONS 

C  ‘  ‘  ‘ 

C 


C 

C 

C@@@@ 

c 


BY  LEAST  SQUARE  METHOD  OF  OBSERVATION  EQUATION 

CALL  INDLSQ  (NC,NR1,NZ1.NUM.KK,K.DIST,TTC0X,TTC0Y.NT.NA,ND.NP, 
NP2 ,  NI , AZ2l  ,AZ34AtfGfUD ,  EdRIDAAftal , STdA ,tflS  , ST0D , 
PR I NTD  PW , Wtf , STDO , FM , BM , BMT , &MTWM , NM , T$ , 

DELTA , NMI ,WKlO , ASNY , ASEX , VM , STDAD) 


WRITE(6,* 
WRITEi  6  * 
WRITE(6  * 


C 

C 


103 


PRINTING  DETAILS 

| S|ANDARD  DEVIATION  OF  UNIT  WEIGHT  =' , STDAD 
'STANDARD  DEVIATION  OF  ADJ.  ANGLES' 

DO  103  1=1. NA 

CALL  RCTmS  (VM(  1.1 ) ,  DUM1D ,  DUM1M , DUM1S) 

DUMMY 1  =  6UMiD*36O0. 0D0  +  DUM1M*60. 0D0  +  DUM1S 
WRITE(6 ,*) 'I  =' , I ,  ADJ.  STDA.  =' , DUMMY 1 

CONTINUE  ' 

WRITE(6,*)' STANDARD  DEVIATION  OF  ADJ  DISTANCES' 
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oooooooooooooo 


c 

c 


DO  104  1=1, ND 
K4  =  NA+I 
WRITE(6,*V 
WRITE(6  *) 
104  CONTINUE 


'I  =  '  ,1,'ADJ.  STDD.  =' , VM(1,K4) 

LUPI  I  ilNUC 

105  FORMAT (I3,1X)I2,1X,I2) 

STOP 
END 


(3  (3  (3(3  (3  (3  (3  (3  (3  (3  (3(30(3(3(3(3(3  (3(3  (3  (3  (3  (3(3  (3  (3  (3  (3(3  (3  (3(3  (3(3(3  (3(9(3 

@@ 

@0 

@0 

SUBROUTINE  LEAST  S( 

1UARE 

@@ 

@@ 

ADJUSTMENT  OF  INDIf 

<ECT 

@@ 

@@ 

OBSERVATIONS 

@@ 

@0 

@0 

@@(3(3  (3  (3  (3(3  (3  (30(3(3  (3  @(3(3(9  (3  (3  (3(3(3  0(3(3  (3(3  (3(3(3  (3  (3  (3  (3(3  (3  (30 


C 

c 


c 

c 

c 

c 


c 

c@@ 

c 


SUBROUTINE  INDLSQ  ( NCI , NR , NZZ , NUMM . KK, K ,  DI ST1 . TCOX . TCOY ,  NT , NA , 

*  ND.tfP.riP2,tfl.AZ2l,Az34.NGRfD.EGGlD,Alta, 

*  stOa.Oisi  !stC)d,pr1ntd.?w,wm  .Stdo.fm, 

*  BM.BHT.BMtWM.NM.TM.DEtTA  NMi.WKld, 

*  ASfsIG  ,A$  EG  .VM;  STAND) 

DOUBLE  PRECISION  DUMMY1 ,NGRID(4J ,EGRID(4) . DUM1D . DUM1M, 

*  DUMlS,XNG(NAT,DfSlLND7.pj'NT>Nt).STD(}CNT), 

*  ASNG^fiPjJ  ASEG(NP)  ,  FM(Nt)  ,BM(^Tt^2)  ,TESTl  ’ 


* 

* 

* 

* 

* 


TEST2.BWTCNP2  NT}:BMTWM(nP2'ntJ,NM(nP2.NP^), 
TM(NP2) .DELTA(NPz) ,NMI(NP2 ,NP2) ’WKlO(Ni) , 
STDA(nA) , STDDCND) '£)UMMY2 ,NUM0,A221 ,AZ34, 
VMfl, NT)  STAND, DI$T1(NA),TC0X(NC1), 

TCOY(NCI) 


I  ,J ,K1, IER,PRINTD, PW. CHECK, NUMM(KK) ,K,KK, 

NMZZ(NA),KON,Ili,Ndl,NNDi 


_ INTEGER 

DATA  Kl/0/,TESTl/0. 00010000000D0/,NUM0/0. ODO/ 


NC4  =  4 

NP5  =  6 

N1  =  1 


20 

30 

C 

C@@ 

C 


SET  WEIGHT  MATRIX 

DO  30  1=1, NT 
DO  20  J=1 ,NT 

WM(I,j)  =  0.  ODO 
IF  (i.EQ.JJ  THEN 
WMil  .JJ  =  1.  ODO 
END  IF 
CONTINUE 
CONTINUE 

FOR  STD.  ANGLE 

IF  (PW.  NE.  1)  THEN 
UO  40  1=1, NA 

DUMMY 1  =  STDA(I) 
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40 


CONTINUE 


CALL  DMSR  (NUM0.NUM0.DUMMY1 .DUMMY2) 

STDO(I)  =  l.ODfl  /  (  (0SIN(DUMMY2)  )**2) 


C@@  FOR  STD.  DISTANCE 
C 

DO  50  1=1, K 

J  =  NA+I  ,  x 
DUMMY 1  =  STDD(I) 

STDO(J)=  l.ODO  /  (DUMMY1**2) 
50  CONTINUE 
C 

COO  SET  UNEQUAL  WEIGHT 

DO  70  1=1, NT 
DO  60  J=1,NT 

"  w  "  THEN 

STDO(I) 


IF  (I.  EQ.  J)  t; 

WM(I,JJ  = 
END  IF 


60  CONTINUE 

70  CONTINUE 
C 

END  IF 
C 

COO  PRINT  WEIGHT  MATRIX 
C 

IF  (PRINTD.NE.  0)  THEN 

WRITE(6,*)  'WEIGHT  MATRIX' 

WRITEC6  *)  '  v 

CALL  USWFM  ( ' R-C. ' ,NC4,WM,NT,NT,NT,NP5) 


C 

C 

< 

C' 

c 


END  IF 

'99  CONTINUE 

K1  =  Kl+1 
C 

IF  (K1  .GT.  2)  THEN 
GO  TO  999 
END  IF 
C 

CHECK  =  0 
C 

C@@  CALL  SUBROUTINE  TO  CALCULATE  "F"  MATRIX 

C 

C 

CALL  CALFM  (NUMM.KK.K.DIST1 .TCOX.TCOY ,AZ21 .AZ34.NA.ND, 

*  NT.NtaRld.teRID.XNG.DfSl.N^ASNd.ASEd.FM) 

C 

c 

c 

c 

C@@  CALL  SUBROUTINE  TO  CALCULATE  "A"  MATRIX  (COEFFICEINT  DX,DY) 

CALL  CALAM  ( NR, NZZ. NUMM.KK.K.DIST1 . TCOX.TCOY, NT .NA.ND.NP, 

*  NP^,NG^ID,EdRld,ASNG,A^EG,BM) 

C 

C@@  TRANSPOSE  OF  "A"  MATRIX 
C 

DO  120  1=1, NT 
DO  110  J=1 , NP2 

c  BMT(J,d)  =  BM(I,J) 

110  CONTINUE 
120  CONTINUE 
C 

COO  "A"  TRANSPOSE  *  W 
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c 

CALL  VMULFF  (BMT,WM,NP2,NT,NT,NP2,NT,BMTWM,NP2, IER) 

"A"  TRANPOSE  *  W  *  A 
C 

CALL  VMULFF  (BMTWM,BM,NP2,NT,NP2,NP2,NT)NM,NP2,IER) 

C3@  "A"  TRANSPOSE  *  W  *  F 
C 

CALL  VMULFF  ( BMTWM , FM , NP2 , NT , N1 , NP2 , NT , TM , NP2 , I ER) 

C 

c 

c 

C@@  (INVERSE  OF  ("A"  TRANSPOSE  W  *  A)  ) 

CALL  LINV2F  (NM.NP2 ,NP2 ,NMI ,N1 ,WK10, IER) 

C 

C@(3  CALCULATE  DELTA  MATRIX  (COORECTION  VECTOR) 

C  (INVERSE  OF  ("A"  TRANSPOSE  W  *  A)  *  "A"  TRANSPOSE  *  W  *  F  ) 
CALL  VMULFF  (NMI ,TM,NP2,NP2 ,N1 ,NP2,NP2 .DELTA ,NP2 , IER) 

C?@  UPDATE  APPROX.  VALUE 

C 

C 

C  KON  =  K-l+NR 


132 

C 

131 

C 

C 


DO  130  1=1, NCI 

Tcoxm  =  Tcoxm  +  delta((i-i)*2  +  n 

TCOY(I)  =  TCOY(I)  +  DELTA((I-1)*2  +  2) 


CONTINUE 


III  =  (K-l)*2 
DO  131  1=1 ,NP 


IF  (NZZ(I)  .EQ.  0)  THEN 

ASEG(I)  =  ASEG(I)  +  DELTA(  (I-l)*2  +  1 
ASNG(I)  =  ASNG(I)  +  DELTA(  (I-l)*2  +  2) 

ELSE 

DO  132  J  =1 , NR 

IF  (NZZ(I)  .EQ.  J)  THEN 

ASEG(x)  =  ASEG(I)  +  DELTA(  II1+(J-1)*2  +1 
ASNGflj  =  ASNG(I)  +  DELTA(  II1+(J-1)*2  +2 
GO  TO  131 
END  IF 
CONTINUE 
END  IF 

CONTINUE 

IF  (PRINTD  .NE.  0)  THEN 
\/RITE(6 , 1020)  K1 

'  :ALL3U$WFM  ( ' R-C.  ' ,NC4 , FM,NT,NT,N1 ,NP5) 

WRI TE( 6 , 1000) 

*RITt  CALL48sWFM  ( ' R-C.  1  ,NC4,BM,NT,NT,NP2,NP5) 

V/RI  TEC  6,10000 

WRITE^cAlL3U$WFM  ('R-C. ' ,NC4)BMT,NP2,NP2,NT,NP5) 
WRITEC6.1000) 

WRI  '  CALL6U3V/FM  (  1  R-C.  1  ,NC4 , BMTWM ,NP2 ,NP2 , NT, NP5) 
WRITEC6.1000) 

WRITE^CALL78sWFM  ( 'R-C.  '  ,NC4 ,NM,NP2 ,NP2 ,NP2 ,NP5) 


J 
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WRITE  6,1000 
WRIT£(6  1080}  ,, 

CALL  USWFM  ('R-C. '  ,NC4,TM,NP2,NP2,N1,NP5) 
WRITE(6 ,1000) 

WRITEC  6  1090) 

CALL  USWFM  ( ' R-C. 1 ,NC4,NMI,NP2,NP2,NP2,NP5) 
WRITE(6 , 1000) 

WRITE(6  1100) 

CALL  USWFM  (R-C. 1 ,NC4, DELTA, NP2,NP2,N1,NP5) 
WRITE(6 , 1000) 

WRITE(6 ; 1110) 

DO  140  1=1 , NP 


140 


CONTINUE 
END  IF 


WRITE(6*,*)  I,'  N  =  ' ,  ASNG(I) 

WRITE(6’*)  f  ' 


E  =  ' ,ASEG( I) 


C 
C 

C@@  CHECK  DELTA  MATRIX 

C 

C 

DO  180  1=1, NP2 

TEST2  =  DABS  (DELTA(I)) 

IF  (  TEST2  .GE.  TEST1)  THEN 
CHECK  =  1 
END  IF 

180  CONTINUE 
C 
C 

IF  (CHECK  .NE.  0)  THEN 
GO  TO  99 
END  IF 
C 

COO  CALCULATE  "A"  *  X  MATRIX 
C 

CALL  VMULFF  (BM,DELTA,NT,NP2,N1 ,NT,NP2,$TD0,NT, IER) 

C@@  CALCULATE  "V"  MATRIX 
C 

DO  190  1=1. NT 

STDO  jf)  =  STDO(I)  -  FM(I) 


190  CONTINU 


C 

C 

C 

C 


139 

149 


VM^l , I )  =  STDO(i; 


CORRECT  OBSERVED  ANGLES  AND  DISTANCES 

DO  139  I  =  1,NA 

ANG(I)  =  ANG(I)+(STDO(I)) 

CONTINUE 
DO  149  1  =  1 .  ND 
NND1  =  NA+1 

DIS1(I)  =  D I S 1 ( I )+STDO(NNDl) 

CONTINUE 

IF  (PRINTD  . NE. 0)  THEN 
WRITE(6 , 1000T 
WRITEC6 . 1120) 

CALL  USWFM  ( 1 R-C. ' ,NC4,STD0,NT,NT,N1,NP5) 
WRITE(6 , 1000)  ' 

WRITEC6 . 1130) 

CALL  USWFM  ( 1 R-C. 1 ,NC4,VM,N1,N1,NT,NP5) 

END  IF 
C 
C 

C@@  CALCULATE  "V"  TRANSPOSE  *  W  *  V  MATRIX 
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CALL  VMULFF  (WM,STD0,NT,NT.N1,NT.NT, FM.NT. IER] 
CALL  VMtH-FF  ( VM, FM ,Nl , Nt ,Ni  ,Nl , Nt , StANO ,Ni ,  IER) 


C 

C 

C 

C 


IF  (PRINTD  .  NE.  0)  THEN 
WRITE(6 , 1000) 

WRITE(6  1140) 

CALL  USWFM  ( 1 R-C. 1 ,NC4 , FM,NT,NT ,N1 ,NP5) 
WRITE( 6 , 1150) 

WRITE(6  1160)  STAND 
END  IF 

CALL  VMULFF  (NMI , BMT.NP2 ,NP2 ,NT ,NP2 ,NP2 , BMTWM.NP2 , IER) 

STAND  =  DSQRT  (  STAND  /  (DFL0AT(NT-NP2))  ) 

DO  240  1=1, NP2 
DO  230  J=1  ,NP2 

NMI  (I,J)  =  STAND  *  (DSQRT(DABS  (NMI(I.J))  )  ) 


END  I 

'230  CONTINUE 
240  CONTINUE 


NMI ( I , J) 


C 
C 

C@@  CALCULATE  "A"*( INVERSE  "A"TRANSPOSE*W*A)  *  "A"TRANSPOSE 

CALL  VMULFF  (BM,BMTWM,NT,NP2,NT,NT)NP2,WM,NT) IER) 

IF  (PRINTD  . NE. 0)  THEN 
wRITE(6, 1000  T 
WRITEC6  1170) 

CALL  US&FM  (‘R-C. 1 ,NC4,WM, NT.NT, NT.NP5) 


C 

C 


END  IF 


DO  280  1=1 , NT 
DO  270  J=1 .NT 

WM  (I,J)  =  STAND  *  (  DSQRT (  DABS(WM( I , J) )  )  ) 


(I.EQ.J)  THEN 
VMi  1, 1 J  = 


IF  ^ _ 

Ml 
END  IF 


WM( I , J) 


270 

280 

C 

999 
C 

C 

1000 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 


CONTINUE 

CONTINUE 


RETURN 


FORMA' 

FORMA- 

FORMA' 

FORMA- 

FORMA' 

FORMA' 

FORMA- 

FORMA' 

FORMA- 

FORMA' 

FORMA- 

FORMA' 

FORMA 

FORMA 

FORMA' 


'W1sl 


1 

/// 


,5 
5X 
5X 
5X 
5X 
5X 
5X 
,  5X 
5X 
5X 
5X 
5X 


7//  5X 
:///,5X 


1  A  MATS T Y 1  /} 

'TRANSPOSE’OF  A  MATRIX' ,/) 

'A  TRANSPOSE  *  W  MATRIX'. 7) 

"  TftlX  , 


a 


'A  TRANSPOSE  *  W  *  A  MATF.._. 

'A  TRANSPOSE  *  W  *  F  MATRIX' 

'THE  INVERSION  OF  A  TRANSPOSE 
'X  MATRIX  OR  DELTA  MATRIX',/) 
'UP-DATE  THE  APPROX.  VALUES',/) 
'V  MATRIX1,/)  ,  x 

'V  TRANSPOSE  MATRIX',/) 

'W  *  V  MATRIX' ,/) 

'V  TRANSPOSE  *  W  *  V  MATRIX' ,/) 


W  *  A' ,/) 
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oooooooooooooo 


1170  mW-Mmto&E  A  TRANSPOSE  *W*A)*A  TRANSPOSE1  ,/) 


END 


(3(9(9|9(3(9(9(B(3(9(9(9(3|3P(3(9P(3(3(9(3(3(9P(3(3@PP(3(3@PP@@PP 
(3(9(9  (9(9(9 

000  SUBROUTINE  FOR  COMPUTING  000 

POO  F  MATRIX  000 

(9(9(9  000 

000000000000000000000000000000000000000 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 


c 

c 


SUBROUTINE  CALFM  (K0UNT,KK,K.DIS,C0X,C0Y,AZ21 .AZ34, 

*  NNA ,  Nfto ,  MnY  ,  GRV ,  GR)(  ,0BANG ,  060IS , 

*  NOS  ASY  ASX  F) 


DOUBLE  PRECISION  GRYjf  4)  .GRX(4)  .OBANGjfNNAl  .OBDIS(NND) , 

N36u , n360R^DUMi ,NUMu , CDl3( 100), TEMPI , 


I NTEGER  1 1 , 1 2 , K , KK , KKK , N , STN , KOUNT( KK) 

DATA  N360/360.  ODO/ , NUMO/O.  ODO/ , N/0/ , STN/1/ 

CALL  DMSR  (N360,NUM0,NUM0,A360R) 

DO  500  1=1, NNA 
DISTAN  =  DIS(I) 

IF  (STN  .  EQ.  1)  THEN 
TEMPI  =  GRXY2) 

TEMP2  =  GRY(2) 

AZB  =  AZ21 
TEMP3  =  ASX(I) 

TEMP4  =  ASY(I) 

CALL  GRIDAZ  (TEMPI ,TEMP2 ,TEMP3 , TEMP4 ,AZF, DUM1 ) 


CDIS(I)  =  DUM1 
IF  (AZF  .GT.  AZB)  THEN 

F( I )  =  OBANG(T  )-(AZF-AZB) 

END  FF1  ^  =  0BANG(I)"(A360R+AZF'AZB) 


N  =  N+l 
GO  TO  400 
END  IF 

IF  (STN  .EQ.  2)  THEN 
TEMPI  =  COX(STN-l) 

TEMP2  =  COY  C  STN— 1 j 
TEMP3  =  GRX( 2) 

TEMP4  =  GRY(2) 

CALL  GRIDAZ  (TEMPI ,TEMP2 ,TEMP3 ,TEMP4, AZB, DUM1) 
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ooooo 


c 

c 


c 

c 


c 

c 


c 

c 


c 

c 


c 

c 


IF 


(KK  .  EQ.  3)  THEN 
TEMP 3  =  GRX '  " 


TEMP4  =  GRY(3) 

END  IF 

TEMP3  =  ASX(I) 

TEMP4  =  ASY(I) 

CALL  GRIDAZ  (TEMPI ,TEMP2,TEMP3 ,TEMP4,AZF,DUM1) 


TO =-Dum 


ELSE 


\ZF.GT.  AZB) 
(I)  =  OBANGC 


AZB)  THEN 
OBANG(I )-(AZF-AZB) 


F(_I)  =  OBANG(I)-(A360R+AZF-AZB) 
END  IF 


N  =  N+l 
GO  TO  400 


END  IF 


IF  (STN  .EQ.  KK)  THEN 
AZF  =  AZ34 
TEMPI  =  GRX( 3) 

'3 


TEMP2  =  GRY 


IF  (DISTAN  .  EQ.  0.  )  THEN 
TEMP3  =  AS X ( I ) 

TEMP4  =  ASY(I) 

ELSE 


TEMP3  =  COX(K-l) 

TEMP4  =  COY(K-l) 

END  IF 

CALL  GRIDAZ  (TEMPI ,TEMP2 ,TEMP3 ,TEMP4 , AZB, DIS1) 


IF  (AZF  .GT.  AZB)  THEN 

T(I)  =  OBANGCI)-(AZF-AZB) 

END  FI(FI)  =  OBANG(I)"(A360R+AZF"AZB) 


N  =  N+l 
GO  TO  400 
END  IF 


IF  (STN  .EQ.  K)  THEN 
TEMPI  =  COX(K-i; 
TEMP2  =  COY(K-l 
TEMP3  =  COX(K-2 
TEMP4  =  COY( K~2 


CALL  GRIDAZ  (TEMP1,TEMP2,TEMP3)TEMP4,AZB,DUM1) 

IF  (DISTAN  .EQ.  0.)  THEN 
TEMP3  =  ASX(I) 

TEMP4  =  ASY(I) 

ELSE 

TEMP3  =  GRX(3) 

TEMP4  =  GRY(3) 

END  IF 

CALL  GRIDAZ  (TEMPI  ,TEMP2 ,TEMP3 ,TEMP4,AZF,DUM1) 
CDIS(I)  =  DUM1 
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c 

c 


c 

c 


c 

c 


c 

c 


400 


C 

C 


c 

c 


IF  (AZF  .GT.  AZB)  THEN 

F(I)  =  OBANGCI)-(AZF-AZB) 

ELSE 

F (I)  =  0BANG(I)-(A360R+AZF-AZB) 


N  =  N+l 
GO  TO  400 


END 

IF 

IF 

(STN 

GT.  2) 

THEN 

TEM 

PI 

=  COX( 

;STN-1) 

TEM 

P2 

=  COY 

STN— 1 1 

TEM 

P3 

=  COX 

'STN-2) 

TEM 

P4 

=  COY( 

;stn-2) 

CALL  i 

GRIDAZ 

(TEMPI 

TEMP3  =  ASX ( I) 
TEMP4  =  ASY(I) 


CALL  GRIDAZ  (TEMPI ,TEMP2 ,TEMP3,TEMP4 , AZF, DUM1) 

CDIS(I)  =  DUM1 
IF  (AZF  .GT.  AZB)  THEN 

F(I)  =  OBANG(T)-(AZF-AZB) 

ELSE 

END  =  OBANG(I)"(A360R+AZF'AZB) 

N  =  N  +  1 


END  IF 


IF  (  (DISTAN  .NE.  0.  )  .AND. 
ADISISTN)  =  CDIS(n 
WRITE(6 ,  )  STN=‘ ,STN , 


END  IF 


KKK  =  KOUNT(STN) 

IF  (N  .LT.  KKK)  THEN 
STN  =  STN 
ELSE 

STN  =  STN+1 
N  =0 
END  IF 


(STN  .  LE.  K)  )  THEN 
CDI S= ' ,CDIS( I ) 


500  CONTINUE 


II  =  NNA 
DO  600  12=1, NND 
II  =  11+1 

Ff  II)  =  OBDI S( 1 2)  -  ADI S ( 1 2 ) 
600  CONTINUE 

RETURN 

END 
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(3(B5?(3!2|?!?!c!pP!S(3P(3(3P(5(?!3(3(3(3(3(3(3(3(3(3!3(3(3!3!3(3(3(3(3(3 


(p@ 

59  SUBROUTINE  FOR  CALCULATING  (3? 

(3?  A  MATRIX  (3? 

(3?  (3(3 


£3  ?  (3  (3  (3(3  @(3(3  (3  (3  (9(350(3(30(3  (3(3(3  00  (3(<3(3(3  (3  (3  (3  (3(3(3  (3(3(3  (3  [3 


SUBROUTINE  CALAM  ( N R , NZ , COUNT , KK , K , DI SS  , COXX , COYY , NT , NA , 
ND,NP,NP2,KrtY,lM,ANY;ANX,AM) 


DOUBLE  PRECISION  KNY(4) ,KNX(4)ANY(NP),ANX(NP_),  DISTAN, 
AM^NT  ,TjP2)iADI^S(  NA)  ,  CdXX(K) ,  CQYY(  K) , 


DUl  ,DU2,Du3,DU4,DUS,DU6 


INTEGER 


1. 11,12,13,14  15  J1.J2, K, KK.K1.K2.K3, 

K4  ,N ,  NN ,  STfj  ,T^M  ,f  EMl  ,  T^m2  ,  J j.jjl ,  jJ2 . 
NOM , COUNT(K^) , Tt , TT1  TEM3 , NU$ED , NZ( NA) , 
NR, L, 110 


DATA  STN/1/ , N/l/ , NN/0/ , TT/0/ , J J/l/ 


NP1  =  NP2-1 

K1  =  NA-1 

K2  =  K-i 

K3  =  K+l 

110  =  ( ( K-l ) *2)  -  1 

(3  (3  (3  (3(3  (3  (3  (5  (3(3(3  (3  (3  (3  (3  (3(3  (3(3  (3(3  (3  (3  (3(3  (3  (3  (3(3(3 

(3 . (3 

(3  COEFFICEINT  OF  COORDINATES  9 
(3  FOR  ANGLE  CONDITION  (3 

( 3  (3 

(3(3(3  (3(3(3!3(3  (3(3(355(3(3  (3(3  (3(3(3(3  (3  [3  (3(3(3  (3  (3(3(3(3 


c 

c 

c 

c 

c 


231 

C 


C 

C 

232 


DO  800  I  =  1 , NA 

DISTAN  =  DISS(I) 

DO  700  J1  =  1 ,NP1 ,2 
J2  =  Jl+1 
AM(I,J1)  =  0.0D0 
AMC I  * J2)  =  0.  0D0 

IF  (STN  . EQ.  1)  THEN 


IF  (DISTAN  .  EQ.  0)  THEN 

DO  231  L=1 ,NR 

IF  (NZ(  I )  .  EQ,  L)  THEN 
JJ  =  I 10+2*L 
GO  TO  232 
END  IF 
CONTINUE 

ELSE 

JJ  =  1 
i EM  =  JJ 
GO  TO  232 
END  IF 


IF  (J1  .  EQ.  JJ)  THEN 


ooo  o  o  o  o  o  oo  o  o  ooo 


c 

c 


DU1  =  ANY 
DU2  =  ANX1 . 
DU3  =  (DUl**2) 


2) 


AM 

AM 


H; 


Jl)  =  DU1/DU3 


J2 


-DU2/DU3 


N  =  N+l 
GO  TO  700 
END  IF 

GO  TO  700 
END  IF 

******************** 


IF  (STN  .  EQ.  2)  THEN 


IF 


(Jl 

OU1 


_J)  THEN 

r(I)  -  COYY ( STN— 1 
((I)  -  COXXfSTN 


DU2 

DU3 

DU4 

DU5 

DU6 


. EQ.  JJ 
=  ANY' 

=  ANX 

=  ^NY(2)  -  C0YY(STH-1 

2  - 


•-U 


2)+(DU2**2) 

-  COYY(STN-l) 
=  KNX(2)  -  COXX(STN-l) 
=  (DU4**2)+(DU5^*2) 


AM 

AM 


(I,J1)  =  -(DU1/DU3)+(DU4/DU6) 
( I , J2;  =  (DU2/DU3)+C  DU5/DU6) 


233 


IF  (DISTAN  .EQ.  0)  THEN 
OO  233  L  =  LNR 


"IF  (  NZ(f)  .  EQ.  L)  THEN 
JJ1  =  I10+2*L 


GO  TO  700 
END  IF 


CONTINUE 
ELSE 

JJ1  =  TEM  + 
TEM1  =  JJ1 
GO  TO  700 
END  IF 
GO  TO  700 
END  IF 


IF  (Jl  .  EQ. 
AM( I ,J1 

“ir’~ 


JJl)  THEN 
=  (DU1/DU3) 


AM(i;j2j  =-fbu2/DU3) 


N  =  N+l 
GO  TO  700 
END  IF 

GO  TO  700 
END  IF 

**************** 


IF  (  (STN  .GT.  2)  .AND.  (STN  .  LE.  K2)  ) 
K4  =  STN  -  2 


IF  (Jl  .  EQ.  JJ 
DU1  =  ANY 
DU2  =  ANX 
DU3  =  (DU 


THEN 

I)  -  COYY(STN-l) 
I)  -  COXX(STN-l) 
t*2)+(DU2*>‘"" 


DU4  =  C0YY(K4)  -COYY(Itn-I) 
DU5  =  C0XX(K4)  -  COXX(STN-l) 


THEN 


78 


nno 


c 

c 

c 

c 

c 


234 

C 


c 

c 

c 

c 

c 

c 


c 

c 

c 


DU6  =  (DU4**2)+(DU5**2) 


AM 

AM 


(I,J1)  =  -DU4/DU6 
(I,J2)  =  DU5/DU6 


JJ1  =  TEM1 
GO  TO  700 
END  IF 

IF  (J1  .EQ.  JJ1)  THEN 


AM 

AM 


(I,J1)  =  -(DU1/DU3) 
(I,J2)  =  (DU2/DU3) 


(DU4/DU6) 

(DU5/DU6) 


IF  (DISTAN  .EQ.  0)  THEN 
DO  234  L  =  1 , NR 

IF  ( N Z ( I )  . EQ.  L)  THEN 
JJ2  =  T10+2*L 
GO  TO  700 
END  IF 

CONTINUE 

ELSE 

JJ2  =  TEM1  +  2 
TEM2  =  JJ2 
GO  TO  700 
END  IF 


GO  TO  700 
END  IF 

IF  (J1  .EQ.  JJ2)  THEN 

DU1/DU3 


c  I  ,Jl)  = 

( I ,  J2)  =  -DU2/DU3 


N  =  N+l 
GO  TO  700 
END  IF 
GO  TO  700 


END  IF 


******************* 


IF  (STN  . EQ.  K)  THEN 
IF  (J1  .  EQ,  JJ)  THEN 

DUl  =  COYYfK-2)  -  COYY(K-l) 
DU2  =  COXX(K-2)  -  COYYfK-1) 
DU3  =  (DU1**2)  +  (DU2**2) 


AM 

AM 


(I,J1)  =  -DU1/DU3 
(I,J2)  =  DU2/DU3 


JJ1  =  TEM1 
GO  TO  700 
END  IF 

IF  (J1  .  EQ,  JJ1)  THEN 

IF  (DISTAN  .  EQ.  0.  )  THEN 
DU4  =  ANY(I)  -  COYY( 

anx(i)  -  COXXf 
"  J4*^2)  +  (DU 


ELSE 


DU5  = 
DU6  = 


TEMP2  =  J 


DU4  = 
DU5  = 
DU6  = 


K-l) 
K-l} 
5**2) 


-  COYY( 

-  COXXf 
2)  +  (DU 


K-l) 

K-l) 

5**2) 
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000 


c 

c 


c 

c 


235 


C 

C 

C 

C. 

C 

c 

c 

c 

c 


700 


AM( I  ,J1)  =  -(DU4/DU6)  +  (DU1/DU3) 
AM(I,J2)  =  (DU5/DU6)  -  (DU2/DU3) 

N  =  N  +  1 

NUSED  =  J1 

JJ2  =  0 
GO  TO  700 
END  IF 
NN  =  NN  +  1 


AM( I , Jl)  =  -(DU4/DU6)  +  ( 
‘  ,J2)  =  (DU5/DU6)  -  ( 


AM(  I 


DU1/DU3) 

DU2/DU3) 


DO  235  L  =  1 .NR 

IF  (  NZ( f)  .  EQ,  L)  THEN 
JJ2  =  I10+2*L 
GO  TO  700 
END  IF 
CONTINUE 


GO  TO  700 
END  IF 

IF  (Jl  .EQ.  JJ2)  THEN 
AM(I,J1)  =  DU4/DU6 
AM( I , J2;  =  -DU5/DU6 

N  =  N  +  1 
GO  TO  700 
END  IF 
GO  TO  700 


END  IF 


****************** 


IF  (STN  .  EQ,  KK)  THEN 
IF  (Jl  .  EQ,  JJj  THEN 

I)U1  =  COYYrK-1)  -  

DU2  =  COXX(K-l)  -  KNX(3l 
DU3  =  (DU1**2)  +  (DU2**^) 


-  KNY(3' 

<  : 


AM 

AM 


(I,J1)  =  -(DU1/DU3) 
(I!J2)  =  (DU2/DU3) 


N  =  N  +  1 
GO  TO  700 
END  IF 

GO  TO  700 
END  IF 

CONTINUE 

NOM  =  COUNT(STN) 

IF  (N  .  LE.  NOM)  THEN 
STN=  STN 
ELSE 

STN  =  STN  +  1 
JJ  =  TEM 
N  =  1 
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ooooooo 


GT.  3)  THEN 
=  TEMI 


IF  (STN 
~  JJ 

TEMI  =  TEM2 
END  IF 
END  IF 


800  CONTINUE 


(?  COEFF.  OF  COORDINATES  0 
@  FOR  DISTANCE  CONDITION  @ 

(3(3(3(3(3(3(B(B(3(3(3(B(3(9(3(B(3(B(9(B(3(3(3P(9P 


C 

C 


c 

c 


c 

c 

c 


DO  3000  II  =  l.ND 

12  =  NA  +  fl 

13  =  II  -  1 

14  =  13  +  (13  -  1) 

15  =  13  +  (13  +  1) 

DO  2000  J1  =  1.NP1.2 

J2  =  J1  +  1 
AM(I2,J1)  =  0. 0D0 
AM(I2  J2)  =  0.  0D0 


IF  (II 
IF 


.EQ.  1) 

VuVii 


THEN 
l)  THEN 

1NYC2)  -  COYY(I) 
DU2  =  KNXf2)  -  COXXfl) 
DU3  =  DSQRT  (  DU1**2  + 

AM( 12 , Jl)  =  — ( DU2/DU3 ) 
AM(  12  ’  v32 )  =  -(DU1/DU3) 


DU2**2) 


END  IF 
GO  TO  2000 
END  IF 


IF  (II 
IF 


.EQ.  ND)  THEN 
(Jl  . EQ.  NUSED)  THEN 
DU1  =  COYY(K-I)  -  KNY(3) 
DU2  =  C0XX(K-1]  -  KNX(3) 
DU3  =  DSQRT  (  DU1**2  +  D 
AM(I2,J1)  =  (DU2/DU3) 

AM(I2  J2)  =  (DU1/DU3) 


U2**2) 


END  IF 
GO  TO  2000 
END  IF 


IF 


(14  .EQ,  Jl)  THEN 
DUl  =  COYYCI3)  - 


)  -  COYY(Il) 
)  -  COXX(Il) 
DU1**2  +  DU 


DU2  =  COXXfl 3 
DU3  =  DSQRT  ( 

AM( 12 , Jl )  =  (DU2/DU3) 
AM( 12 ( J2 )  =  (DU1/DU3) 

GO  TO  2000 
END  IF 

IF  (15  .EQ,  Jl)  THEN 
OU1  =  COYYCI3)  -  COYY 
DU2  =  COXX(.I3)  -  COXX 
DU3  =  DSQRT  (  DU1**2 


2**2  ) 


hi) 

+  DU2*y 


2  ) 


i. 
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c 

c 


2000 


-  AM( 12 , Jl)  = 
AM( 12 | J2)  = 

END  IF 

CONTINUE 


-(DU2/DU3) 

-(DU1/DU3) 


C 

C 


3000  CONTINUE 
RETURN 

END 


(3  (3  (3  (3  (3  (3  (3  (3(3  (3(3  (3  (3(3  (3(30(3  (3  (3  (3  (3  (30(3(3  (3  (3  (3  (3@(3(3  (3  (3  (3(3(3@ 

00  00 

(3(3  SUBROUTINE  TO  CALCULATE  GRID  (3(3 

00  AZIMUTH  AND  DISTANCE  (3(3 

(3(3  (3(3 

000000000000000000000000000000000000000 

SUBROUTINE  GRIDAZ  (XI ,Y1 ,X2,Y2 , ANGLR , DIST12 ) 


REAL*8  X1.X2.Y1 , Y2.ANGLR.ANG0 .ANG90 . ANG180 .ANG270 
REAL*8  AN690ft,A^ll86R,AN2t0R,DiFX  ,DlfY,DISTl2 

DATA  ANG0/0.  0D0/.ANG90/90. 0D0/ .ANG180/180.  0D0/ 
DATA  ANG270/270. 600/ 


ANG90 ,ANG0 ,ANG0 .ANG90R  ) 
ANG186 , ANGd , ANG6 , AN 180R  ) 
ANG270, ANGO  ANGO, AN270R  ) 


CALL  DMSR 
CALL  DMSR 
CALL  DMSR 

DIFX  =  X2-X1 
DIFY  =  Y2-Y1 

IF((DIFX.  EQ.ANG0).  AND.  (DIFY.  EQ.ANG0))  THEN 
ANGLR  =  ANGO 

ELSE  IFfDIFX. EO.  ANGO)  THEN 
IF(DIFY.GT.  ANGO)  THEN 
ANGLR  =  ANGO 
END  IF 

ANGLR  =  AN180R 


ELSE 


END  IF 


IF(DIFX.GT.ANGO)  THEN 

ANGLR  =  ANG90R-DATAN(DIFY/DIFX) 

ANGLR  =  AN270R-DATAN(DIFY/DIFX) 
END  IF 


DIST12  =  DSQRT(DIFX**2+DIFY**2) 

RETURN 

END 


(3(3(3(3(30000000000000000000000000000000000 
(3(3  (3(3 

00  SUBROUTINE  FOR  CHANGING  THE  (3(3 

@0  ANGLE  FROM  DEGREE,  MINUTE,  AND  (3(3 

(9(3  SECONDS  TO  RADIAN  (3(3 
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00  00 
00000000(300000000000000000000000000(3(300 

SUBROUTINE  DMSR  (  D1,M1,S1,R1  ) 

REAL*8  D1,M1,S1,R1,PI1 , N60 , N3600 , N180 

DATA  N60/60. 0D0/ ,N3600/3600. 0D0/ ,N180/180. 0D0/ 

PI1  =  4. 0D0*DATAN(1. 0D0) 

R1  =  ((D1+(M1/N60)+(S1/N3600))*PI1)/N180 

RETURN 

END 


000000000000000000000000000000000000000 
00  00 

00  SUBROUTINE  FOR  COMPUTING  GRID  00 

00  COORDINATES  FROM  KNOWN  AZIMUTH  00 

00  DISTANCE  AND  STARTING  C00RDI-  00 

00  NATES  00 

00  00 

000000000000000000000000000000000000000 


SUBROUTINE  UTM  (XI ,Y1 ,DIS12,ANG12R,X2,Y2  ) 

REAL*8  X1,Y1,X2,Y2,DIS12 ,ANG12R,DIFX ,DIFY 

DIFX  =  DIS12*(DSIN(ANG12R)) 

DIFY  =  DIS12*(DC0S(ANG12R)) 

X2  =  Xl+DIFX 
Y2  =  Yl+DIFY 

RETURN 

END 


C 

c 

c 


c 


000000000000000000000000000000000000000 
00  00 

00  SUBROUTINE  FOR  CHANGING  THE  00 

00  ANGLE  FROM  RADIANS  TO  DEGREE,  00 

00  MINUTES,  AND  SECONDS  00 

00  00 

000000000000000000000000000000000000000 

SUBROUTINE  RDMS  (  RR1,D2,M2,S2  ) 

REAL*8  RR1,D2,M2,S2,PI1,N60,N180,TDEG,DIF,TMIN 

DATA  N60/60. 0D0/ ,N180/180. 0D0/ 

PI1  =  4. 0D0*DATAN( 1.  0D0) 

TDEG=  (N180*RR1)/Pll 
D2  =  DFLOAT(IDINT(TDEG)) 

DIF  =  TDEG-D2 
TMIN=  DIF*N60 
M2  =  DFLOAT(IDINT(TMIN)) 

DIF  =  TMIN-M2 
S2  =  DIF*N60 

RETURN 

END 
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